Fano Resonances in Quantum Dots: A Non-Perturbative Role for Potential-Like 

Scattering 
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We consider the physics of transport through quantum dots in the presence of two tunneling paths. 
The first path sees electrons hopping on and off the dot while the second path is modeled through 
a potential scattering-like term. To study the effects of potential scattering, we employ a modified 
version of the Anderson model. Such a model can be exactly solved through the Bethe ansatz, thus 
allowing a comprehensive and exact analysis of the zero temperature linear response conductance. 
We find transport properties to be extremely sensitive to the introduction of a potential scattering 
term. Indeed the presence of such a scattering term, inter alia, induces a series of first order quantum 
phase transitions. Focusing on the Kondo regime of the quantum dot, the non-perturbative effect 
of potential-like scattering can be directly tied to both the breaking of particle-hole symmetry and 
the interlinking of charge and spin degrees of freedom in the Anderson model. The sensitivity to 
potential scattering is also reflected in a set of complementary exact diagonalization computations. 
The consequences of this analysis extends to observations in general of Fano resonances in quantum 
dots as well as to the physics of transport through quantum dots embedded in Aharonov-Bohm 
rings. 

PACS numbers: 73.23.Hk, 85.35.Gv, 72.15.Qm, 71.10.Pm 



I. INTRODUCTION 

Fano resonances are a commonplace phenomena aris- 
ing from a competition between two competing scatter- 
ing paths, one path sensitive to the energy of the scat- 
tering state, one path not. These resonances were first 
studied by the eponymous U. Fanoi in the context of 
photoionization experiments in atomic systems where a 
continuum of states interacting with a single resonant 
level serve as the active ionization pathways. This worki 
predicted that interference between two scattering paths 
will in general lead to asymmetrically shaped resonances 
as a function of scattering energy. Research on Fano 
resonances in atomic systems remains a topic of active 
interest Resonances of similar origin have long been 
known to characterize scattering off nuclei^ More re- 
cently they have been observed in the context of inter- 
ference in quantum well systems^ polymer films j£ tun- 
neling into and between multi-wall carbon nanotubcs, 8 
and STM measurements of adsorbed magnetic atoms on 
metallic substratesj2*i2*ii particularly interesting from a 
mesoscopic viewpoint in that Kondo physics comes into 
play. 

In this article we will be interested in yet another man- 
ifestation of Fano resonances, that involving single elec- 
tron transistors (SETs). Interference between two tun- 
neling paths, one path involving an SET embedded in 
an arm of an Aharonov-Bohm ring, was first studied 
in Ref. [H However asymmetric resonances typical of 
the phenomena were only first reported in the conduc- 
tance profiles of SETs several years thereafterii^^ In 
this work, the exact geometry of the interfering tunneling 
paths through the SET was left unspecified, a situation 
again seen in the experiments of Ref. Il6i More recently, 



Kobayashi et aiiZilfL have reported Fano resonances in 
SETs embedded in Aharonov-Bohm rings where neces- 
sarily the two tunneling paths are clearly delineated. 

An SET (or more colloquially, a quantum dot) is fabri- 
cated from a gated semiconductor hctcrostructurc such as 
GaAs/AlGaAs and is formed from a segregated nanome- 
ter scale sized region of a two dimensional electron gas 
(2DEG). Such quantum dots are highly tunable. Both 
the hopping between the dot and the bulk of the 2DEG 
as well as the dot's energy levels relative to the 2DEG's 
Fermi surface can be adjusted through metallic gates 
placed on the hctcrostructurc. At sufficiently small tem- 
peratures (much smaller than the dot level spacing and 
tunneling widths) only the electron levels on the dot near- 
est in energy to the Fermi energy of the bulk 2DEG con- 
tribute to transport, turning the system into a quantum 
impurity problem. As is the case in such problems, the 
effective dimensionality of the bulk 2DEG is d e g = 1 in 
regards to questions concerning the impurity, a result of 
focusing upon low energy impurity S-wave scattering^ 
With d e ff = 1, the bulk electrons can be treated as living 
in leads. 

The remarkable ability to control the various parame- 
ters of a quantum dot has allowed Kondo physics to be 
observed in such a settingi22*2ii22i Tuning the dot such 
that one electron sits on the dot level nearest the Fermi 
surface gives rise to the much of the same physics as 
an isolated magnetic impurity^. With the introduc- 
tion of a second tunneling path, Kondo physics does not 
disappear^ but as we will argue, takes on a considerably 
different guise. 

In Figure 1 we sketch out the relevant energy scales of 
the dot with only its energy level nearest the bulk Fermi 
surface drawn. The energy levels of the two 'leads' are 
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FIG. 1: Sketched are the relevant energy scales of the single 
active dot level. 



marked as forming a continuum filled to levels marked 
by two chemical potentials, [ij, and [Ib.- The active dot 
level is characterized by two energy scales, one marking 
the energy of single occupancy, e^, the other marking the 
energy of double occupancy, U + 2e^. U here represents 
roughly the charging energy incurred from the addition 
of the second electron into the level. Separating the ac- 
tive level from the two leads are two tunnel barriers. As 
drawn, the hybridization of the electrons living on the 
level with those in the leads gives rise to the broadening 
of the level to a width, r. Represented in the capacitively 
coupled gate voltage, V g , is the ability to shift relative 
to /Lt£ and 

The properties of a single level SET are believed to be 
well captured by the Anderson model Theoretical 
computations of the finite temperature linear response 
conductance in the Kondo regime, a non-trivial quantity, 
both using numerical renormalization group techniques^ 
and the Anderson's model underlying integrability^ well 
match the observations m& Moreover qualitative predic- 
tions resulting from the out-of-equilibrium Anderson 
modelj2£ namely the splitting of the zero bias anomaly 
into two peaks, have been observedSiSi (although the 
predicted zero temperature features of the peakaS^^ 
have not yet been seen2&). 

A continuum Anderson model governing the quantum 
dot is given by 

TiAnd. = ^ \ ~i dx(c' a(J (x)d x c a(7 (x)) 

a = L.R 

+V(cl a d cr + dlc aa )\ x=0 [ + e d ^ n <y + Un^n^. (1.1) 

Here c aa /d are electrons living in the leads/dot. Hopping 
between the two is governed by V . As indicated in Figure 
1, e<2 marks the dot single level energy while U describes 
the additional cost of double occupancy. 



This description of an SET yields a single, resonant 
tunneling path for electrons traveling from the right to 
the left lead. A second non-resonant tunneling path in 
the vicinity of the dot can arise in a variety of manners. 
This path can occur through the dot itself^*i& a product 
perhaps of the detailed geometry of the confining region 
of the 2DEGS Alternatively, the second path can be 
spatially separated as occurs when a quantum dot is em- 
bedded in an arm of an Aharonov-Bohm ringfiSiiLiS as 
pictured in Figure 2. This is in some sense preferable as 
the two tunneling paths are then marked concretely by 
the two arms of the ring. 

We model the second tunneling path in the most mini- 
mal way possible by adding to the Anderson model above 
the term, 



SH = V L r {c La c RlT + h.c. 



(1.2) 



x=0 



Such a term has been considered previously in a number 
of places^2^i*22iS It allows the electrons to bypass the 
dot through a direct lead-lead tunneling term. In adding 
(57i, we allow the possibility of multiple windings as elec- 
trons transit successively from one lead to the other via 
different paths. This modified Hamiltonian is thus ap- 
propriate to describing an Aharonov-Bohm ring in a two 
terminal geometry^SiiLiS not the four terminal setup of 
Ref. |3 The second tunneling path is potential-like in 
structure. Indeed with the introduction of even and odd 
combinations of electrons, 
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the addition is precisely a potential scattering term. 

At the centre of this work lies the discovery that 
this Hamiltonian, Hau<\. + 8H- i« integrablc. While it 
has long been known that the Anderson model alone 
is integrablej2£i2& only recently has it been realized that 
the addition, SH, to the Hamiltonian does not spoil the 
integrability^ As will be shown here, the perturbation, 
87i, modifies the integrablc structure of the Anderson 
model in a straightforward way. While this modification 
is minor in a technical sense, it will have profound effects 
on the transport properties of the system. 

At the heart of the model's integrability lies the knowl- 
edge of the model's exact eigenfunctions. From the 
knowledge of the exact eigenfunctions, we are able to 
extract the behaviour of excitations above an interact- 
ing Fermi sea. In particular we can compute the dressed 
scattering matrices of the excitations both with them- 
selves and off the impurity. While the integrability of 
the Anderson model has been long established, it was 
only recently understood^ how to wed integrable data 
with a Landaucr-Buttiker framework to access transport 
properties. In particular, in Ref. |2^ it was shown how 
the knowledge of the excitations above the Fermi sea and 
their associated scattering phases allow the computation 
of both the in-equilibrium and out-of-equilibrium con- 
ductances. 
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FIG. 2: Drawn is a ring connected to two leads. In the lower 
arm of the ring a quantum dot is embedded. Electrons transit 
from lead L to lead R by either hopping on and off the dot 
embedded in the lower arm of the ring or by traveling through 
the upper arm of the ring. Electrons may also move from lead 
L to lead R by executing multiple windings about the ring. 



Much of the content of Ref. |25| was devoted to under- 
standing how to treat both finite temperature, T, and 
voltage, V, transport in the context of integrability. At 
finite T and V various uncertainties plagued the discus- 
sion. Away from the Fermi surface, electronic excita- 
tions could be constructed only in an ambiguous fash- 
ion. The association of a finite energy electron with 
a particular integrable excitation was then left unfixed. 
This ambiguity was overcome only in a number of cases 
with some educated guesswork. In this way Ref. |25| 
was able to reproduce the finite temperature linear re- 
sponse conductance determined from numerical renor- 
malization group calculations^ Here, in contrast, we will 
focus upon zero temperature, linear response transport. 
There is no ambiguity in identifying the necessary elec- 
tronic excitations and their scattering phases at the zero- 
temperature Fermi surface. As such the results contained 
herein should be considered exact and without approxi- 
mation. 

At zero temperature, the linear response conductance 
in the Vlr — Anderson model, 7iAnd., is described 
by the Friedel sum rule. The Fricdel sum rule relates 
the scattering phase, S e y^, at the Fermi surface to the 
number of electrons sitting on the dot, 



= 2im 



dVl- 



(1.3) 



With the second tunneling path activated, we find that 
a variant of the Friedel sum rule still holds. But now 
the scattering at the Fermi surface is related to the total 
number of electrons displaced by the impurity (here the 
impurity should be thought of as both the dot and 5TL): 



27TO; 



impT/i 



2tan~ 1 (V p ). 



where V p = Vlr/2. While the additive constant, 
2tan _1 (VJ,), is important, the crucial difference between 
the two expressions for the scattering phase lies in that 
rii m p is distinct from rid- At least for U > 0, exact diago- 
nalization computations that complement the integrable 
analysis indicate that n. 



impT/i 



involves both contributions 



from the dot degrees of freedom as well as bulk electrons 
in the leads: 



^imp^ 



Mmp4 



dx 



dx 



(cl^(x)c e l(x)} - pbulk 
(cl L (x)c el (x)) - pbulk 



(1.5) 



Here pbuik equals the unperturbed (absent the impurity) 
density of states of the lead electrons. 

It is hardly surprising that deviations of the electron 
density in the leads should affect scattering at the Fermi 
surface. This was specifically forewarned as a possibil- 
ity in the proof of the Fricdel sum rule for the Anderson 
model by LangrcthiiS And we know such deviations in 
the lead electron density arc important for determining 
the finite field scattering phase (and so the magnetocon- 
ductance) in the Kondo modeli^iS Instead what turns 
out here to be surprising is that nj mp = ni mp f + ^i m p| is 
not a continuous function of V p . 

As a function of V p , nj mp sees two discontinuities. The 
first discontinuity occurs at V p non-zero, i.e. at V p = 
0*. The origin of this discontinuity is perhaps clearest at 
6d = —U /2, at what would be the particle-hole symmetric 
point of the Anderson model if V p = 0. V p finite then 
breaks particle-hole symmetry. In this breaking of the 
symmetry, a two-fold degeneracy of the ground state is 
broken leading to a first order quantum phase transition. 
The two putative ground states carry different electron 
occupancy. Their breaking leads then to a discontinuous 
change in nj mp by one electron. At least at the point 
6d = — U/2, this conclusion is supported by an exact 
diagonalization study. 

A second first order transition function of V p oc- 
curs approximately when the two tunneling paths become 
equal in strength. The sign of V p at which the transition 
occurs is a function of 6d- If 6d > — U/2, transport is 
naturally expressed in terms of particles (as opposed to 
holes) and the transition occurs at V p = V p crlt > where 
y crit is given by (for r < U) 



crit 



2r 
77' 



(1.6) 



This positive value of V p marks the energy necessary to 
ionize an electron from the vicinity of the impurity. If 
6d < —U/2, transport is most naturally given in terms of 



holes and the second transition occurs at V v 



< 



(1.4) 



'p ^ -p 

0. V p = — V p rlt then marks the energy necessary to ionize 
a hole from the impurity site. 

These two transitions in V p do not persist to all e^. 
Rather they are bounded (roughly) to fall in the range 
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FIG. 3: A phase diagram of the Anderson model as a function of ed and V p . The horizontal dashed (blue) lines mark first order 
quantum phase transitions while the vertical solid (red) lines mark both the termini of the phase transitions as a function of 
td and regions of rapid variation in nimp- 



-U 



<e d <e; l 



where e™' 



equals 



By making the identifications, 



crit 



£ cnt mar ]j S the point where the disturbance in the impu- 
rity density of states due to V p occurs above the Fermi 
surface and so is excluded from contributing to the scat- 
tering phase. 

These transitions are encoded in Figure 3, a phase di- 
agram of the Anderson model as a function of ed and V p . 
Each area of the phase diagram is marked by its charac- 
teristic value of rii m p. The horizontal dashed (blue) line 
segments mark the two sets of first order phase transi- 
tion. On the phase transition lines themselves, rii mp ~ 1 
(not marked in Figure 3). The termini of these phase 
transitions are delineated by solid vertical (red) lines. Al- 
though the phase transitions at V p = and V p = V^ crlt 
are drawn so as to end both simultaneously at e d rlt and 
— U — e" 1 ', the termini differ by a value of order T. The 
vertical lines also mark regions of rapid (although not 
discontinuous) variation in 7ii mp . 

The conductance, G, is simply expressible in terms of 
the scattering phase (at H = where S e -\ = S e i = 5 e ), 



G 



2^sin»(i(f. 



2tan- 1 (V p )). 



(1.8) 



e = cot(i(<5 e + 2tan- 1 (F p ))); 

q = -cot(2tan -1 (Vp)), 
G can be put into a Fano-like form, 



G 



(e + q) 



r,Y2 



h 1 



Vp 2 e 2 



1 



(1.9) 



(1.10) 



However this form is somewhat misleading, e is not be 
identified with e^/r as is true in the U = case. In 
fact unlike the non-interacting case, we find the many- 
body correlations introduced by finite U can lead to the 
non-vanishing of G as ed is varied with q still real. It 
would be interesting to reconsider then the analysis of 
the Fano resonances found in Ref. (although the large 
magnetic fields at which the resonances were observed 
might largely eliminate many-body effects). 

The origin of the phase transitions and their atten- 
dant consequences can be understood not only as a con- 
sequence of many-body interactions but of mere quantum 
mechanics as well. The exact solution of the V p ^ An- 
derson model dictates that the ground state consists of a 
sea of two-electron repulsively interacting bound states. 
(This feature of the ground state is independent of V p 
and holds true for the V D = Anderson model.) While 
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this feature results from the exact treatment of the many 
body physics in the problem, the origin of the discontin- 
uous behaviour in ni mp can be seen in the bare bound 
states themselves. If one computes the finite V p scat- 
tering phase of a two-electron bound state off the im- 
purity (a quantum mechanical problem), one sees that 
it is characterized as a function of energy by a number 
of non-analyticitics. The close relationship between the 
undressed bound state scattering phase, the impurity oc- 
cupancy, n lmp , and so ultimately the dressed scattering 
phase at the Fermi surface results in transport proper- 
ties sharing these same non-analyticities. The explicit 
construction of the bound state wave function and a dis- 
play of the non-analyticities are carried out in Section 
II.B.3. 

The appearance of non-analyticities in the bound state 
scattering phase and so transport properties is in part a 
consequence of the treatment of the two tunneling paths, 

Htun = V ^ ( C L^ + dl c l<r) U=0 

l = L,R 
<r=t,l 

+w P ( c L c ^ + 4x c ^)u=o- (i-ii) 

In Wtun, the two tunneling terms are presented as both 
infinitely narrow with delta-function strength. How- 
ever in any real experimental system the two tunnel- 
ing terms will have finite, differing widths. In general 
unequal widths can lead to changes in the physics. If 
for example the width of the lead-lead term is much 
greater than the width of the dot-lead coupling, i.e. 
wicad-icad ~> Wdot-ioad, the effect of the second tunnel- 
ing path on the physics becomes relatively trivial. The 
results presented in this article are presumed valid only 
if ifiead-icad < ^Mot-lead- More is said on this in Section 
II.B.I. 

Htun also treats the two tunnelings as occurring at the 
same point, x = 0. If in the effective theory we were to 
displace the scattering terms by some arbitrarily small 
distance the non-analyticities disappear. The distinction 
can be made clear by considering the single electron U = 
scattering phase. The scattering phase of c e = (cl + 
cr)/V2 at wavevector q is given by 

<5 e (<?) =-2tan- 1 (— !— +V P ). (1.12) 

q - £d 

Notice that V p appears within the argument of the arc- 
tan. In the case of a quantum dot embedded in an 
Aharonov-Bohm ring, this form of the scattering phase 
results from topology, from the possibility of multiple 
windings around the ring. While we arc ultimately inter- 
ested in the scattering phase of bound states, this form 
of the scattering phase is crucial if non-analyticities in V p 
are to arise. Now if instead we shift (say) the potential 
scattering to x = e, i.e. 

L = L. R 



+ 2V p Y ( c L c ^+ c L c ^)u=c (i.i3) 

the scattering phase reads instead 

5 e (q) = -2tan _1 (^— ) - 2 tan -1 (V p ) . (1.14) 

q - td 

We now have excluded the possibility of multiple wind- 
ings and the limit of the bound state scattering phase as 
a function of V p will be entirely well behaved. 

Generically the tunneling terms will not exactly coin- 
cide and we will have the situation as given in i|1.14f> . 
This is reassuring. If potential scattering was always to 
have drastic consequences we would need to explain why 
its presence is not always felt in experiments in mag- 
netic impurities in bulk metals. The effective Anderson 
Hamiltonian of a magnetic impurity in a bulk metal will 
generally not possess particle-hole symmetry. Thus it 
will include, however small, a potential scattering term. 
However the coincidence of dot and potential scatter- 
ing centers as given in requires some enforcement 
mechanism (such as perhaps the topological possibility 
of multiple windings arising for a quantum dot embed- 
ded in a ring). This enforcement mechanism will not be 
present in the general case and potential scattering will 
be innocuous. 

What then is the likelihood that the scattering phase 
appearing in (|1.12|) is relevant to understanding Fano res- 
onances in quantum dots? Certainly this phase is generic 
theoreticallyi22i2i The differences between our work and 
other analyses detailed in this article arise not because of 
the form of S e (q), but because of the treatment of inter- 
actions on the quantum dot. But what is the situation 
experimentally? While Fano resonances have been ob- 
served, their observation does not necessarily preclude 
the form of the scattering phase in Eqn. I|1.14fl . O nc 
possible means of distinguishing between the two forms 
of the scattering phase lies in the effect of finite inter- 
actions. As we present in Section VI, various features 
of the observations found in Ref. ^| suggest a role for 
finite interactions only consistent with 11.12f) . However 
this cannot be considered definitive while the role inter- 
actions play in other reports of Fano resonancesiSiiLi& is 
undetermined. 

The dramatic effect potential scattering presents is not 
only a consequence of the coincidence of scattering cen- 
ters. It reflects the entanglement of charge and spin de- 
grees of freedom that occurs in the Anderson model. In 
the Kondo model, one may add a potential scattering 
term with strength, V p , and maintain the model's inte- 
grability. If the scattering phase at the Fermi surface is 
then computed onc finds 

s B m = * + f(Vp)> (i-i5) 

where n is the original contribution of the screened spin 
to the scattering phase while / is a smooth function of 
the strength, V p . This scattering phase then behaves in 
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the same fashion as the scattering phase for an Anderson 
model with potential scattering but without coincident 
scattering centers. The trivial effect of potential scat- 
tering is a result of the complete uncoupling of spin and 
charge degrees of freedom in the Kondo model. Poten- 
tial scattering acts directly only in the charge sector while 
all the non-pcrturbative physics occurs in the spin sector. 
Thus potential scattering is unable to effect the Kondo 
physics in an interesting fashion. This limitation is re- 
moved however for the Anderson model where the two 
sectors can interact. 

A consequence of this behaviour is that the Ander- 
son model with potential scattering cannot be treated 
via a Schrieffer- Wolff (SW) transformation^ To imple- 
ment the standard SW transformation on 7iAnd. + STi, 
one should first transform away the potential scattering 
term via 

c e (x ^ 0) = e-' lS / 2 6(-x)£ e (x) + e u/2 d(x)c e (x); 
c e (0) = cos{5/2)5(0) + VV p cos 2 (5/2)d; 

5 = -2tan _1 (t^). (1.16) 

This gauge transformation transforms Ti. e into the ordi- 
nary Anderson model with renormalizcd parameters, id 
and V: 

id = e d - F 2 F P cos 2 (5/2); 

V = Vcos(ty2). (1.17) 

The SW transformation taking the Anderson model to a 
Kondo model goes forward as usual by freezing out charge 
excitations. The scattering of the transformed electrons 
at the Fermi surface is the Kondo scattering phase 

S = #Kondo = 1", (1-18) 

while the scattering phase of the original electrons will 
have the form 

S = S Koado -2tB,n- x (V p ), (1.19) 

that is, one restores the scattering phase due to potential 
scattering that was originally transformed away. Com- 
paring this result with (|1 . 14|> or (|1.15|) . suggests that an 
SW transformation performed in such a fashion implic- 
itly assumes the potential scattering term is not exactly 
at x = 0. Alternatively, the gauge transformation in 
(I1.16|) together with the SW transformation seem to mis- 
treat the coupling between the spin and charge sectors of 
the theory. We will sharpen our critique of the gauge 
transformation in (|1.16|) in what follows (see in particu- 
lar Section V). 

While the SW transformation fails to describe the 
physics of the V p ^ Anderson model, this failure does 
not imply that Kondo physics is absent. One hallmark 
of Kondo physics is the presence of the Abrikosov-Suhl 
resonance. This is an enhancement in the low energy 



spectral weight of the impurity degrees of freedom. This 
ensures quantities such as the zero temperature magneto- 
conductance vary rapidly on the order of the Kondo tem- 
perature, Xfc. At V p 0, low energy spectral weight is in 
general still present in the Kondo regime (e<2 ~ —[7/2) of 
the Anderson model but in modified form. With V p finite, 
spectral weight present near the Fermi surface is shifted. 
The nature of this shift can be read off from the value 
of riimp in Figure 3. For V p = we have n lmv ~ 1 for 
£d ~ —[7/2. Taking V p finite, nj mp shifts to ~ or 2. A 
shift of nj mp to 2 coincides with low lying spectral weight 
just above the Fermi surface being transferred to just be- 
low. The corresponding absence of spectral weight at low 
but positive energies leads to the magnetoconductancc 
being insensitive to small changes in H. Analogously, a 
shift of riimp from 1 to implies that low energy weight 
has been transferred from just below the Fermi surface 
to just above. Here in contrast the magnetoconductancc 
becomes all the more sensitive to small changes in H. 
The magnetoconductance is discussed in detail in Sec- 
tion II. B. 

With these comments in hand, we now outline the 
course of the paper. A preliminary, considerably more 
terse, treatment of this subject has already appeared^! 
In Section II we outline the exact solution of the finite 
V p ^ Anderson model. We begin by introducing the 
model followed by the development of the Bethe ansatz 
solution in several steps. By constructing the single and 
two-particle eigenstates of the model wc are able to ex- 
tract the scattering matrices/phases. The structure of 
these scattering matrices demonstrate that the model is 
integrable. Beyond the scattering in the system, the two 
particle eigenstates allow us to exhibit the nature of the 
bound states. As we have already indicated, much of the 
physics is encoded in the scattering of the bound states 
off the impurity. By further sketching the development of 
multi-particle eigenstates, we are able to argue that in all 
details but the impurity scattering phase, the equations 
describing the quantization of momentum in the system, 
the Bethe ansatz equations, are identical to the V p = 
case. This allows us readily to develop the method by 
which dressed scattering phases of electrons at the Fermi 
surface are computed. In the first part of Section II wc 
limit ourselves to the case td > —[7/2. We so end Section 
II by extending the analysis to td < —[7/2. This exten- 
sion will prove essential when we come to compute the 
conductance exactly at td = —[7/2. 

In Section III we develop the resulting properties of 
the linear response conductance. Wc begin by comput- 
ing the zero field conductance. We carefully associate 
discontinuous changes in rii mp with similar changes in 
the conductance. The appearance of the discontinuous 
changes finds support in exact diagonalization compu- 
tations which show similar discontinuous changes in the 
ground state occupancy. Key features of G as both a 
function of V p and td are described here. In particular 
the ingredients going into the phase diagram of Figure 
3 are established. Having described G(H = 0), we then 
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FIG. 4: The unfolding of the left (L)/right (R) leads into 
chiral fermions extending from x = — oo to x = oo. 



turn to understanding the magnetoconductance. Here 
we are able to characterize analytically the behaviour 
of the H 7^ conductance at = — U/2. Away from 
6d = — U/2 we rely upon numerical solutions of the equa- 
tions describing the dressed scattering phases. 

The role two-electron bound states play in the exact 
solution of the V p ^ Anderson model is ultimately a 
consequence of a finite U. With U = 0, the behaviour 
of the conductance. G, as a function of V p is relatively 
bland. This points to a non-pcrturbative role for U. To 
reinforce this, in Section IV we consider the solution of 
the V p ^ Anderson model for U < 0. With U < 0, all 
non-analyticities in V p disappear suggesting that various 
impurity quantities do not possess well-behaved power 
series expansions in U/T as with the V p = model. 

In the penultimate Section V, we analyze an alterna- 
tive approach to this problem, an approach predicated on 
the use of Dyson equations and developed in Ref. |3(] and 
Ref. 001 This approach leads to transport quantities be- 
ing analytic in both U and V p and so to considerably dif- 
ferent results than presented in this article. We will argue 
that the differences arise because the Dyson equations 
treat the impurity occupancy, 7ij mp , as equivalent to rid, 
the number of electrons on the dot itself. We will equate 
the Dyson equations to the results obtained through ap- 
plying the gauge transformation in i|l.lbf> . This will en- 
able us to argue that even at the quantum mechanical 
level, the Dyson equations seem to incorrectly treat finite 
interactions. We believe that the origin of the discrep- 
ancy between the two approaches lies in that the Dyson 
equations presume implicitly both V p and U are pertur- 
bative quantities and so end with quantities well behaved 
in V p and U. 

In Section VI, the final section, we consider the obser- 
vations reported in Gores et al>i^ as these served as the 
initial motivation for this work. While hardly determina- 
tive, some generic features present in these experiments 
are well explained using the formalism developed in this 
work. 



II. EXACT SOLUTION 

A. Description of System 

We use a more general Hamiltonian than in and 
(|1.2|> to describe the lead-impurity system, 

"H = H i ca ds + #tun + ^tun + ^dot ! 

/oo 
dx ( ~ ic\ a (x)d x ci a (x))\ 
-oo 

H Ln = ^(4<A + 4<v)|x=o; 

l = L,R 
<r=T4 

H tun = Vlr (c{ a c R(J + c R(J c La )\ x=0 ; 

<T=t,l 

H lei = 2J (Vllc^clv + V RR c\ ia c R(7 )\ x= o; 
H do t = £d ^2 n a + Un^ii. (2.1) 

«r=T4 

Here again cj '/ci„ are electron operators in the leads 
while d^/dcr are the corresponding operators on the dot 
with n a = S a d a . The leads in Figure 4 are semi- infinite 
and have electrons which move both to the left and the 
right. By "unfolding" the leads (in the continuum limit), 
each lead can be thought to run from x = — ootoa; = oo 
but with chiral electrons which are solely right moving. 

The path from L to R taking the electron through the 
dot is modeled by H^ un . We have allowed for the possi- 
bility that tunneling between the dot and the two leads 
is asymmetric, i.e. Vl ^ V R . The second path from L 
to R (for example, the upper arm of the ring in Figure 
2) is modeled by H^ un , a potential scattering- like term. 
We have also included a term in the Hamiltonian, i? re f, 
by which electrons again experience a single-body poten- 
tial but one which does not scatter the electrons into the 
other lead. Although likely to be generically present, we 
here consider such a term for technical reasons alone, as 
will be seen shortly. 

Hdot describes the energies associated with the dot. 
The dot single particle Hilbert space is two fold: one for 
occupancy by a spin ^-electron and one for occupancy by 
a spin- J, electron. The chemical potential of these levels 
is controlled by n a while Coulomb interactions on 

the dot are represented by the term Un^n^. This is in 
fact the only interacting term in the problem - the leads 
arc treated as Fermi liquids. This is reasonable as the 
one-dimensional leads ultimately arise out of a dimen- 
sional reduction of electrons living in two or three di- 
mensions. Although U is the only non-trivial interaction 
in the problem, it presents the same set of complications 
seen in a fully interacting system. 

To treat this Hamiltonian we want first to map it to a 
model involving a single electron specie. This cannot be 
done for the general values of Vll, VrR: Vlr, and Vl/Vr. 



8 



However if we restrict ourselves to the parameter space, 

V LR 



V ll = -Vrr 



2V l Vr 



(2.2) 



we can proceed. Note that if the dot-lead coupling is sym- 
metric (or if Vlr = 0) then Vll and Vrr must vanish. 
With this constraint on the parameter space we introduce 
even and odd electrons, c e / Q , via 



c e /o = (Vl/rCl ± Vr /l cr)/^Jv[+V~1 
The Hamiltonian in this basis then reads 
H = H e + H , 

where 



(2.3) 



(2.4) 



H P = 



+{2T) 1 l 2 (c\ a d a + 4c ea )\ x=0 + 2V p c\ a c ea \ x=0 
+e d ^ n <y + Un- ] n i ; 



<T=1,l 



H = ^2 J dx - icl a (x)d x c oa -(x) - 2V p cl a c or7 \x=o- 



(2.5) 



The parameters T and V p are given by 



2 

TV L r 
2V l V r 



(2.6) 



r represents the total broadening of the electron dot level 
due to coupling to the leads while V p represent the effec- 
tive potential scattering seen by the even and odd elec- 
trons. The odd sector of the theory, uncoupled to the 
dot, is trivial. The even sector, we will argue, is solvable 
via Bethe ansatz. 

Having undertaken this recasting of the Hamiltonian 
we must describe how to maintain contact with the orig- 
inal problem: the transport of electrons from the left 
lead to the right lead and vice versa. We can do this 
as follows. Let Tw be the transmission amplitude of an 
electron from lead I to lead I' and TZi the reflection am- 
plitude of an electron in lead I back into the same lead. 
In the even and odd basis it will turn out we can read- 
ily compute the scattering phases, S e /S 0l of an even/odd 
electron off the impurity (here the scattering must be a 
pure phase as the even/odd sectors are uncoupled and 
both chiral). T/1Z can be readily expressed in terms of 



V R e 



V L e l 



VrTZr + VlTlr; 
VlTIr - VrTlr. 



(2.7) 



Similar relations hold with Trl and TZl- Thus 



1~lr 
TIr 



V l Vr 



VI 



vi 



(e 



iS r . 



y2 e iS e + y2 e U 



vi 



vi 



(2.8) 



By parity, we must have \Tlr\ = \Trl\ and \1Zr\ = \1Zl\- 
The zero temperature linear response conductance is then 
given by 

p 2 p 2 V 2 V 2 ft — ft 

G = 2 e -\T LR f = 2^^sin 2 (^) (2.9) 

For most of the paper we will assume Vl — Vr. Only 
when compa ring theory with the experimental observa- 
tions of Ref. [lj will we relax this constraint. 

The prescription at finite energy/finite temperature for 
relating S e /S to T and 1Z is in general fraught with am- 
biguities. This is discussed in scholastic detail in Ref. I2M 
In this paper we will focus upon T = linear response 
quantities. Taking this limit lifts these ambiguities mak- 
ing the above transformation precise. 



B. Construction of Multiparticle Eigenstates for TL e 

1. Single Particle Eigenstates 

We first examine the single particle eigenstates of H e ■ 
This will give the bare scattering phase off the impurity. 
We want to look for states of the form 



IVv) 



dx{g a (x)c i lJ {x)} + e„d\ 



|0>. (2.10) 



Solving Schrodinger's equation, Helper) = <7|Vv), yields 
the following constraints upon g a {x) and e a : 



9. 



(7 (0)(2r) 1 / 2 + e (ieo 



~id x g a (x) + 8(x)((2r) 1 / 2 e a + 2V p g a (0)) - qg a {x). 

(2.11) 

Taking g a {x) to be of the form 

g a (x) = 6{x)e lqx+l5 / 2 + 6(-x)e lqx ~ lS/2 ; 

g a (0) = cos(%)/2), (2.12) 
we can readily solve (I2.11|) finding expressions for e a and 



(2T) 



1/2 cos(%)/2) 



q - £d 

5(q) = -2tan- 1 (—^— + V p 

q- e<z 



(2.13) 



Note that the potential scattering strength appears inside 
the argument of the tan" 1 . This, ultimately, is a reflec- 
tion of the coupling between the spin and charge sectors 
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of the Anderson model (and if appropriate, the possibil- 
ity of multiple windings about an Aharonov-Bohm ring) . 
It thus sets the scene for the non-trivial physics of the 
potential scattering term. 

S(q) contains information on the non-interacting (U = 
0) states of the retarded dot Greens function: 

ImGrH = ^-d q S(q)\ q=u 

Z7T 



( q - ~e d y + i* 



(2.14) 



where T and id are the renormalized values of the dot 
level broadening and dot chemical potential in the pres- 
ence of a finite V„ 



(d - V p f. (2.15) 

We see the presence of the potential scattering term 
serves to narrow the dot resonance as well as shifting 
the position of the energy level. In terms of correlators 
involving dot degrees of freedom alone, these renormal- 
izations, we believe, will be the sole response to the po- 
tential scattering. However the correlators relevant to 
transport go beyond dot correlations. In these correla- 
tors the effects of potential scattering will be non-trivial. 

The resulting single particle phase in (|2.13|) results 
from treating the lead-lead coupling, V p , as well as the 
dot-lead coupling, T, as delta- function scatterers. How- 
ever as pointed out in the introduction, in any exper- 
imental system the two couplings will have some finite 
widths, Wdot-icad and u>i e ad-iead- It is relatively e asy to 
show that in order to sec the physics encoded in (|2.13|) 
we require u>dot-icad > wicad-icad- To see an indication 
of this, suppose we broaden out the lead-lead scattering 
term, replacing V p 8(x) in (|2.11|) with V p w(x) where w(x) 
is some even smooth function with unit weight. The re- 
sulting solution to the modified (|2.11|) is 



9a (x) 



exp i 



dx'2V p w{x')) 



x {e{x)e lqx+lS/2 + 9(-x] 

c -a/„ (2r) i/2 Cos(<W2) 
q — td 

-2tan^ 1 (- 



Jqx-i8/2 



q - £d 

Thus the effective scattering phase of an electron as it 
traverses from x = — oo to x = oo is 

S cS (q) = -2tan -1 (- 1 



\-2V p 

q- Q 

Unlike with the electron scattering phase in H2.13[) . we 
would expect all finite U transport properties to be 
smooth functions of V p . Then in order to see the effects 
predicated upon i|2.13[l . we (at least) require w do t-ioad > 

Wload-lead- 



2. Two-Particle Eigenstates 

The examination of the two-particle eigenstates will al- 
low us to construct the S-matrix for the model. In prin- 
ciple we will then be ready to compute the multi-particle 
eigenstates. The only non-trivial two-particle states are 
those with total S z = 0. We thus focus on these. Such 
states have the form 



m = 



dxidx 2 g(xi , x 2 )ct {xi)c\ (x 2 ) 



dxe(x)(c\(x)d\ - c\{x)d\) + fd\d\ 



|0). (2.16) 



Solving the Schrodinger equation, TL e \^) — E\ip), leaves 
us with the following set of consistency equations, 

Eg(xi,x 2 ) = -i(d Xl + d X2 )g(xi,x 2 ) 

+ v / 2f(<S(x2)e(a;i) + <5(x 1 )e(x 2 )) 

+ V p g(xi, x 2 )(S(xi) + 6(x 2 )); 

Ee(x) = -id x e(x)+V2Tg(x,0)+V2TfS(x) 

+ e d e(x) + 8(x)V p e(0); 

Ef = 2V2T e(0) + 2e d f + Uf, (2.17) 

together with assuming g(xi,x 2 ) is symmetric. To solve 
these relations, we employ the following ansatz 

3(^1^2) = g q {xi)g p {x 2 )(j){xi - x 2 ) 

+g P (xi)g q {x 2 )(j>{x 2 - xi); 

e(x) = e P 4>(x)g q (x) + e q cj>(-x)g p {x); 

E = q+p. (2.18) 

Here g q / p (x) and e q / p are the coefficients appearing in 
the one particle wavefunction corresponding to energies 
q/p (i.e. (|2.12() and 12.1311 ). cj)(x) is an expression of 
the phase change incurred when two electrons are inter- 
changed, x\ <-> x 2 . We assume it takes the form 



4>(x) = 1 + ia(p, q)sign(x). 
Under this ansatz, we find 



(2.19) 



/ = 
a(p,q) = 



-2V2T 
1 



e p cos(S(q)/2) + e q cos(S(p)/2) 



2e d + U- (q + p) 
2UT 



q - p 2e d + U - q - p 



(2.20) 



From a(q,p) we can determine the scattering matrix for 
two electrons. 

This S-matrix, Sv\al(p, q), where o\ and a 2 are the 
spins of ingoing electrons while a' x and a' 2 are the spins 
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FIG. 5: Plots of the bound state impurity scattering phase as 
a function of energy of the bound state. We take U + 2ed = 0, 
U = 1 and T = .1. 



of the outgoing electrons. This S-matrix, on the grounds 
of SU (2) invariance, must have the form 



Salal — b(j>, 



c(p, q)P^ 



(2.21) 



where I, 



P 



(71 (72 

6?' 



1 



5%\8al is the identity matrix while 
is the permutation matrix. The coef- 



ficients, b(p, q) and c(p, q), are related to a(p, q) via 

(xi-x 2 ) l + ia(p,q) 



b(p, q) - c(p, q) = 



4>{x 2 -xi) l-ia(p,q) 

g(p)-g{q) + i 
g(p) - g{q) - i' 



(2.22) 



where g(p) = (p-e d - U/2) 2 /(2UT). With this relation 
in hand, b and c can then be fixed by considering the 
scattering of two electrons with equal spin. As two such 
electrons do not interact, necessarily S%% = 1. Hence 



b(p,q) +c(p,q) = 1. 



(2.23) 



So then 



b(p, q) 

c(p, q) 



g(p) ~ g(g) 
g(p) - g{q) 

— i 



(2.24) 



g(p) - g(q) - i 

What is to be noticed is that this is exactly the same S- 
matrix as for the Anderson model with V p = 0. The sole 
change V p finite introduces to the analysis is a change in 
the impurity scattering phase. 



FIG. 6: Plots of the derivative with respect to energy of the 
bound state impurity scattering phase. We take U + 2ed = 0, 
U = 1 and T = .1. 



3. Two-Particle Boundstates 

In constructing the two-particle eigenstates we deter- 
mined the S-matrix of scattering between opposite spins 
of momentum p and q. In particular the state in (|2.16(l 
is a singlet with an associated scattering phase, 



5 s i 



singlet 



b(p, q) - c(p, q) 



g(p) - g{q) + i 



(2.25) 



gip) - q(q) -i' 

This S-matrix has poles in it corresponding to the choice 

gip) = A + i/2; g(q) = A - i/2, 

where A is for now an arbitrary real parameter. At 
these values, the pole in 5 s i ng iet has a positive imaginary 
residue indicative of a stable bound state. 

For U > 0, the ground state is constructed out of such 
bound states. Much of the physics is thus contained in 
the behaviour of the bound state scattering from the im- 
purity. As such it is worthwhile to construct the bound 
state wave function explicitly. We assume the bound 
state has the same form as l|2.16ll . The various compo- 
nents of the eigenstate, g(xi, x 2 ), e(x), and /, still satisfy 
(|2.17|) . However instead of (|2.18|l . we take as an ansatz 
for these functions the following 

g(x!,x 2 ) = 9(x 1 -x 2 )h(xi,X2)f+(x 1 )f-(x 1 ) 

+ 6(x 2 ~ x 1 )h(x 2 ,xi)f + (x 2 )f-(xi); 
e{x) = e-9(x)h{x,0)f+(x) 

+ e+e{-x)h{0,x)f-{x); (2.26) 



where 



f±(x) = 9{x)e l5 ±l 2 + 6{-x)e- l5 ±' 2 - 
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h(xx,x 2 ) 
S± 

e± 
E 



_ p ix-(xi+X2) — y(xi,X2). 



-2tan _1 ( 

x ± iy - e d 

2Tcos((5±/2)^ 



Vn 



x±iy- e d 



2x. 



(2.27) 



Using e(x = 0) = (e_/ + (0) + e+/_(0))/2, one can readily 
show this ansatz provides a solution to l|2.17[) provided x 
and y, the real and complex parts of the momentum of 
the two particles forming the bound state, satisfy, TU + 
2y(x - U/2 - e d ) = 0. 

The scattering phase of the bound state off the impu- 
rity is given by 



-'bd. state 



= s 



1 bg f ix - e d - iy)(l - iV p ) - ^^ 
i \{x - e d - iy)(l + iV p ) + iT 



X 



(x- e d + iy)(l - iV p ) - iT 
(x-e d + iy)(l + iV p )+iT 



(2.28) 



We plot the bound state scattering phase as function of 
the energy E = 2x. At V p = the scattering phase un- 
dergoes a discontinuous jump at E w —U. This jump is 
of magnitude 7r. It however does not in the end produce 
any discontinuous behaviour as the derivative of 5 

bd. state 

with respect to the energy, E, is continuous. This in turn 
will imply that the bare impurity density of states asso- 



J 



ciated with the bound state is smooth. However small 
changes in V p from introduce radical changes. If V p > 0, 
the size of the discontinuity jumps to 2ir and <9E<5bd. state 
sees a large positive variation in the vicinity of E as —U. 
As V p — * + , this variation becomes a delta function, 
S(E + U). If V p < 0, the Vp = discontinuity in <5bd. state 
is smoothed out as evidenced in Figure 5. However this 
smoothing also leads to a large but negative variation in 
$bd. state (Figure 6). As V p — > 0~, this variation becomes 
approximately — 5{E + U). 

These large variations in the scattering phase of the 
bound states are determinative of the physics, having 
profound effects on the linear response conductance. It 
is worthwhile to stress then that the origin of any novel 
behaviour in the linear response conductance is found at 
the level of quantum mechanics. When we construct the 
ground state out of bound states, this behaviour is only 
modified quantitatively not qualitatively. The dressing 
of the scattering phase due to interactions between the 
bound states does not fundamentally modify the discon- 
tinuity structure of the bound state scattering phase. 



4- Multiparticle Eigenstates 

Multiparticlc eigenstates see the same phenomena: the 
role V p plays is solely encoded in the single particle 
scattering phase. To see this, we write down an N- 
particle eigenstate with M j electrons. We will suppose 
N > M > 2. Similar arguments work for the other cases. 
This wavefunction has the general form 



\iPn,m) = I dxi ■ ■ -dx N g(xi, . . . ,x N )c\(xi) ■ ■ ■ ^(x M )c\{x M +i) • • -c\(x N ) 

dxi ■ ■ -dxM-idxM+i ■ ■ -dx N ei(xi, . . . ,x M -i,x M +i, ■ ■ ■ ,x N ) 

xcj(xi) • ■ ■ c\(x M -i)c\{x M +i) ■ ■ -c\{x N )d\ 



j dx\ ■ ■ ■ dxudxM+i ■ ■ ■ dxN-i e]{x\, . 

fdx 1 ...dx M _ 1 dx M+1 ...dx N - 1 /;,-, 



i XM,XM+1, ■ ■ ■ , XN-l) 

xcj(xi) ■ • • c\(xm)c\(xm+i) ■ ■ ■ c^(x N -i)d\ 

■ , XM-1,XM+1, ■ ■ ■ : %N-l) 

xcj(xi) • ■ ■ c\(x M -i)c\{xm +1 ) ■ ■ ■ c\(x N _i)d{d|. (2.29) 



We will look for eigenstates with total energy E = then take the coefficients g(xi),ei(xi), e^(xi), and f(xt) 
J2iLi 1i where the q t are single particle energies. We will to be of the general form 



g(xi,...,x N ) = ^2 Y\_9q aW (x t ) ■ A N , M (xt\a); 

ctESn 



12 



ei(xi, ... ,xm-x,xm+\, ■ ■ ■ ,xn) 



e-f(a;i, . .., xm,xm+u ■ ■ ■ , x N -\) 



f(xi, . . .,XM-1,XM+1, ■ ■ ■ , %N-l) 



M e <?-<™> II 9q^{xi) ■ B lNM (xi,x M = 0|cr); 

a£S N *=1 

N-l 

(TV — M) ^ e q„ [m) II 9q aW {xi) ■ B rNtM (xi,x N = Ojcr); 

ct^Sn J — 1 

7V-1 

(TV-M)M e 9CT(A/) e Mjv) n^w^O 



Here 2^<t6S n * s a sum over P ermu tations of TV. <7 g (x) 
is the electron portion of the single particle eigenstate 
defined in l|2.10[) with energy q. Similarly e q is the dot 
portion of this same wavefunction. These are the sole 
ingredients of the wavefunction with an explicit depen- 
dence upon V p . A^M^N-M(xi\a) are coefficients which de- 
pend upon both a and the relative ordering of the Xj (i.e. 



r 



i=l 



xC N M (xi,x M =x N = 0|cr) 



(2.30) 



A^m^n-m only changes value as the Xi cross one another). 
Their forms are analogous to those of the Hubbard model 
detailed in Ref. 0- Moreover and most importantly, the 
A^m^n-m have no explicit dependence upon V p and are 
identical in form to those of the V p = Anderson model. 

To see this, we analyze the TV-particle Schrodinger 
equation. In doing so we find pieces of the form 



N 



= 



= 



( - id X( + V p S(x l ) \g(xi, ■ ■ ■ ,x N ) H ; 

i=l ^ ' 

^ (- id Xi + VpS(xi)jei(xi, ■ ■ ■ , x M -i, x M +i, ■ • ■ ,x N ) H ; 

Y y ~ + V p<>( x i)j e t( x l, ■ ■ ■ ,X M ,X M +1, ■ ■ ■ ,X N -l) H ; 

Y ( - i dxi + V p 6(xi))f(xi,---,x M -i,XM+i,---,XN-i) + 



i = l 



(2.31) 



r 



But because we have constructed g, e±, e-f, and / out of 
products of single particle cigenfunctions and g q (x) has 
the property 



d x g q {x) = V p g q (0)6(x) 



(2.32) 



all terms proportional to V p cancel straightforwardly. 
The equations (or pieces thereof) that remain to be 
solved are then equivalent to the ordinary Anderson 
modeliS*^ Although a V p dependence in these equations 
lurks in g q (x) and e q , it never has to be explicitly exhib- 
ited in their solution. 



C. Bethe Ansatz Equations: Analysis for td > —U/2 

In the previous section, we argued TV— particle wave- 
functions of energy E — 53i 9« could be explicitly con- 



structed. Moreover the S-matrix describing the scatter- 
ing of two electrons was shown to be identical to that 
of the V p = Anderson model. V p appears only in the 
impurity scattering phase. This implies that the quanti- 
zation conditions (the Bethe ansatz) for the qi when the 
system is placed in a system of finite length, L, are near 
identical to those of the usual Anderson model. 

On the basis of the above discussion, the Bethe ansatz 
equations for this system, mimicking Ref. l36l take the 
form 



N 

n 



JqjL+i5(qj,V p ) 



^q^g(gj) + V 2 



-tt g(gj) ~ Ac +i/2 
ll ig ( qj )-\ a -i/2> 



M 

n 

/3=1 



A/3 + i . 
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5(q,V p ) = -2tan~ x ( 

q — £d 

(q-e d - U/2f 



9(q) 



2UT 



V 9 ); 
(2.33) 



These equations depend upon V p only through S(q,V p ). 
In fixing the allowed values of qj, they depend upon a 
set of M rapidities, {A Q } equaling the number of spin 
| electrons in the system. The total spin, S z , of the 
wavefunction is 2S Z = N — 2M. 

These equations have a rich structure of solutions for 
{<?j } and {A Q }. But at T = only two types of solutions 
are important: 

i) qj a real number; 
and 

ii) A a a real number to which are associated two 
complex q's, given by 



<?(<£) 
x(X) 



g(x(X a ) =F iy(X a )) = X a ± i/2 
U/2 + e d - VUf(X 

VuT(-x + t/x?~ 



1/4) 1 / 2 . 



1/4) 



1/2. 



These two solutions are the only ones appearing in 
the ground state. x(X) and y(X) defined above satisfy 
TU + 2(x — U/2 — €d)y, the condition identified in the 
construction of the bound state wave function in Section 
II.B.3. 



alter the spin, q is known as a charge excitation and p(q) 
the charge distribution. They are derived by taking log's 
of the Bethe ansatz equations, 



M 



L + S( qj ,V p ) = 2irN J -Y^0i{g(q J )-Xp) ] 



(3=1 



M N-2M 

2TrJ a + ^2d 2 (X a -X l3 ) + ^ ( A « - 9{Qj)) 

,3=1 i=i 

= -2Lx(X a ) - 2ReS(x(X a ) + iy(X a ), V p ); 



Ax) = 2tan~ 1 (-a;) - 7T. (2.34) 
n 



where Nj and J a are quantum numbers with the J a sat- 
isfying 

N-2M N - 2M 
2— <J ^^2— 

By taking appropriate derivatives in (|2.34() 



p{l) = ^ + + 9'(q) dX ai (g(q) A)a(A); 

<t(A) = _ f Q dX'a 2 (X' - X)a(X') 



7T L 

B 

dqai(X - g{q))p(q), 

D 



(2.35) 



D. The Ground State 



where 



The ground state of the finite V p Anderson model con- 
sists of the same set of solutions as at V p = 0. As such a 
ground state with spin 2S Z — N — 2M consists of 

i) N-2M real ^'s; 

ii) M real A Q 's and the 2M associated q a ±- 

The construction of this ground state holds for > 
—U/2. We can understand the case Ed < —U/2 through 
an explicit construction focusing upon holes rather than 
particles. There will be call to do so. We can alterna- 
tively perform a particle-hole transformation. This trans- 
formation however changes the sign of V p . We will find 
reason to examine both approaches. 

In analyzing the ground state we want to consider the 
thermodynamic limit. Thus we no longer treat qj and X a 
as discrete but rather derive smooth distributions (per 
unit length), p(q) for the qj's and cr(A) for the A a 's and 
their associated complex g's. These distributions, p(q) 
and cr(A), mark out how the solutions of the Bethe ansatz 
equations are distributed when the system is in its ground 
state. As each A in the ground state reduces the ground 
state spin by 1 (i.e. 2S Z — N — 2M), a A is known as a 
spin excitation and <r(A) the spin distribution. As each 
q in the ground state at fixed particle number does not 



A(g) = ^d q S(q); 

A (A) = --d x Re6(x(X)+iy(X),V p ); 

IT 

a n {x) = —d x 9 n (x) = — 

2ir it (n z + ix z ) 



(2.36) 



In the equations for the charge and spin distributions 
appear the Fermi surfaces, Q and B, and "band" bottoms 
Q and —D. These mark out the range of q and A over 
which excitations appear in the ground state. —D is the 
lower band edge of the charge excitations. As each A has 
two associated complex q's, we expect q q+ (Q) + q q _ (Q) — 
2x(Q) = —2D, thus determining Q~ The Fermi surfaces 
Q and B arc fixed by insisting that the overall spin and 
particle number, 



N-2M = L dqp{q); 

J-D 

M = L dXa(X), 



(2.37) 



are reproduced. 
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E. Scattering Phase at the Fermi Surface: T = 
Transport 

We are now in a position to access the quantity we are 
ultimately interested in, the zero temperature scattering 
phase, <5 e , of electrons at the Fermi surface. To compute 
8 e we employ an argument used by Andrei (Ref. 14(1) in 
the computation of the magnetoresistance in the Kondo 
model. We imagine adding an electron to the system. Its 
momentum will be quantized according to p = 2im/L. If 
the impurity was absent, the momentum would be deter- 
mined by the bulk properties of the system and we would 
write p = Pbuik- In the presence of the impurity, the to- 
tal momentum of the electron, p, experiences a shift away 
frompbuik ) necessarily scaling as 1/L, 

P = Pbulk +Pimp/L. 



dqa x {X - g{q))pbuik(q), 



(2.40) 



D 



and 



Pimp(g) = A(q)+g'(q) / dX ai(g(q) - A)cr imp (A); 

Jq 

rQ 

01mp(A) = A (A) - / dX a 2 (A - A)cr iinp (A ) 
Jq 

dq ai {X-g{q))p imp {q). (2.41) 

-D 

We now relate the impurity density to the impurity mo- 
menta and so the impurity scattering phases. 

The total momenta of the q and A excitations are deriv- 
able from the continuum limit of logarithms of the Bethe 



5 e is then identified with p imp . In the context of the ansatz equationsfup to possible terms independent of A 
Bethe ansatz, p imp is readily computed. an d q - see Ref. 1251) : 

As we arc interested in quantities scaling as 1/L, it is 
useful to isolate the portions of the densities, p(q) and 
cr(A), scaling as such. With this in mind wc write 



1 

L' 
1 

T 



P{l) = Pbulk(<?) + TPimp(<2); 



<T(A) = Obulk(A) + — (Timp(A). (2.38) 



Here pimp(g) and o"i mp (A) determine the total number 
of electrons displaced by the impurity, rii mp , that is the 
electrons sitting on the dot, n^, plus the deviation of 
electron density in the leads away from their equilibrium 
value of pbuik: 



^impf 



I'd] 



dx 



[cl l (x)c e1 (x)) - p hul] 



dAer imp (A) + I dqp imp (q); 

-D 



^impj. 



n dl + I dx 

Q 



c ll( x ) c el(x)) ~ Pbulk 

dXa imp (X). 



(2.39) 



Typically the deviation of the electron density in the 
leads is taken to be zero. However this is not appro- 
priate in a magnetic field nor will it be appropriate here. 
Indeed this deviation will be of the same magnitude as 
the dot electron number /ndl • 

To derive equations governing the densities, 
Pbuik (g)/P imp(g) and fTbuik(A)/fT imp (A), we substi- 
tute (|2.38() into l|2.35|l and segregate the bulk terms from 
those scaling as 1/L 



Pbuik(g) = ^-+g'(q) / dA ai(g(q) - A)cr bu ik(A); 
27r Jq 

a bul k(A) = -f^l- / dX'a 2 (X' - A)a bulk (A') 



p(q) = q + S(q)+ / dX a(X)9 1 (g(q) ~ X); 

Jq 

p{\) = 2x{X) + 2ReS(x(X) + iy(X)) 

rQ f-B 
+ dX'a(X')6 2 (X-X')+ dqp(q)0 1 (X-g(q)). 
Jq J-d 

(2.42) 

The impurity momenta correspond to the portion of the 
total momenta scaling as 1/L: 

rQ 

Pimp(<7) = <%)+ / dX (T imp (X)6 1 (g(q) - X): 
Jq 

Pimp {X) = 2Re6{x(X)+iy{X)) 
rQ rB 
+ dX'a imp {X')8 2 (X-X')+ t% imp (g)0i(A -<?(<?)). 
Jq J-d 

(2.43) 

Wc note that 

Pimp(g^-oo) = -2tan~ 1 (V p ); 

p imp (A^oc) = ~4tan- 1 (V/ p ). (2.44) 

and so at the band edges, the impurity momenta reduce 
to their bare values. 

The impurity density and momentum are intimately 
related. It is straightforwardly checked that 

d q Pimp{q) = 2iTp(q); d\p imp (X) = ~2tto(X). (2.45) 

These relations will ultimately allow us to express <5 e j n 
in terms of JT-impt/^impi; that is, verify a version of the 
Friedel sum holds in this case. We again stress that in 
this proof of the Friedel sum rule, 5 e ^ / j is not related to 
merely the electron occupancy of the dot, but the total 
electron number displaced by the impurity. 
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In this description of the system an electron is a com- 
posite object. It arises as a gluing together of a charge 
excitation, q, as well as a spin excitation, A. To see this, 
imagine a spin | electron being added to the system. 
The number of electrons in the system as a whole goes 
from N to N + 1 with the system spin increasing by 1/2. 
This means we must add a single q excitation to the sys- 
tem. But in doing so, the range of quantum numbers, 
— N ~2M < j q < N-2M available ^ the X a increases by 
1. This additional slot is in effect a hole in the spin dis- 
tribution. Thus a spin | electron is composed of a charge 
fc-excitation and a spin A-hole. 

The electron scattering phase is the sum of the scat- 
tering phase of these two elementary objects 

Si = Pimp(g) + (Pimp(A) - Pimp (A = Q)). (2.46) 

The additional constant, p(X = Q), ensures that at the 
band edge (q = —D and A = Q) the electron scattering 
phase equals its bare value of — 2 tan _1 (V^). We obvi- 
ously have some degree of freedom in how we choose q 
and A. But we want to construct an excitation at the 
Fermi surface. To do so, we must choose q and A to be 
at their respect Fermi surfaces, q = B and A = Q. Us- 
ing relations l|2.46[) . we can then write 8\ in terms of the 
impurity densities: 



I dqp imp (q) + 2n / d\a imp (X) 



-D) 



27rn impT - 2tan~ 1 (V^). 



(2.47) 



We have thus demonstrated 8\ is (up to an additive con- 
stant) proportional to the number of electrons displaced 
by the impurity. 

To determine the scattering of a spin J, electron we 
employ particle-hole symmetry. A particle-hole transfor- 
mation is implemented via 







-<?); 


- 


■* c t(~ 


<?); 


d\ - 


■* d i; 












-> -U 


- £d 


v P - 


- -v p 





(2.48) 

Thus the scattering phase of a spin J. hole is related to 
that of a spin | electron via 



S l ho (-U-e d ,-V p )=6Ue d ,V p ). 
Thus at the Fermi surface we have 



(2.49) 



SL(-U-e d ,-V p ) = 2tt [ Q d\a imp 

Jo 



(A) 



+2ir dqp imp {q)~ 2 tan" 1 (V p ) 

J-D 

= 27rn im p T (ed) - 2ta n - 1 (U p ). (2.50) 

As n impT (e d , V p ) = -n impl (-U-e d , -V p ) mod 1, we have 
(up to a factor of 2-k) 

SL(- U ~ e d> ~ V p) = -27rn imp j,(-J7 - e d ,-V p ) 

-2tan~ 1 (U p ). (2.51) 

At the Fermi surface 5 l ho (-U - e d ,-V p ) = -S^(-U - 
e d , — V p ). Hence we have shown the Friedel sum rule holds 
for spin | electrons. 

Imagine now that we are at the point where U +2e d = 
and are working at zero field, H = (thus N = 2M and 
B = 0). Hence spin j and | electrons scattering identi- 
cally. Moreover we expect the phase to be insensitive to 
the sign of V p due to particle-hole symmetry. 

Sl(U = -2e d , V p ) mod27r = 5\{U = -2e d , -V p ). (2.52) 

We will however find this to be not naively the case. This 
seeming violation of particle-hole symmetry at the point 
U + 2e d = is a result of the Bethe ansatz yielding a 
state other than the ground state over a range of V p . In 
fact we are able to construct two candidate ground states 
at the point U + 2e d . As we vary V p , the true ground 
state varies between these two possibilities. Once this is 
factored into the calculation, particle-hole symmetry is 
restored. Which ground state to use will be discussed 
in some detail when we discuss the computation of the 
linear response conductance. 

To construct the second candidate ground state, we 
must begin with the Bethe ansatz solution from a dif- 
ferent starting point, one in which we focus on hole not 
particle excitations. We do so now. 



F. Bethe Ansatz Equations: Analysis for < —(7/2 

In the previous sections we analyzed the Bethe ansatz 
equations and ultimately determined the scattering of 
electrons at T = at the Fermi surface. But this anal- 
ysis was limited to the case of U + 2e d > 0. Here we 
consider the case U < ~2e d . Our primary interest in 
doing so is to understand how particle-hole symmetry is 
preserved at the point U — 2e d = 0. To consider U strictly 
smaller than —2e d , we could simply exploit a particle-hole 
transformation. 

To analyze the case U < —2e d , we recast the even 
Hamiltonian in (|2.5|) as 

H e = I J dx - ic ea (x)d x c\ a {x) 

-V2T {d a c\ a + c ea d\ )\ x= q + V p c err cl a \ x= o^ 
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where 



-e d h ° + Unify + V P + 2e d + U, (2.53) 
«r=U 



(d = -f-d - U ; 

Vn = ~Vn. 



(2.54) 



In this form we can construct hole states in exactly the 
same fashion as particle states governed by the parame- 
ters (V p , id)- For example the single hole states take the 
form 



hole 



dx{g a (x)c a (x)} + e a d a 



where here the vacuum state |0) is completely full, i.e. 
cj.|0) = 0. For a single hole of energy q we find 

9a (x) = 6{x)e lqx+l5/2 +9(-x)e lqx - ,s / 2 ; 

r 



%) 



-2tan _1 (- 
-2tan _1 ( 



q — ed 

r 



v P ). 



(2.56) 



The bulk S-matrix describing two-hole scattering is deriv- 
able in a similar fashion to as before and is given by 



^CTi (To 



b(p, q) 

c(p, q) 



Kp> q)lZl% + c (Pi q) p °l°l 

g(p) - g(q) 
g(p) - g(q) -i' 

—i 

g{p) - g{q) - 



(2.57) 



where g(p) = (p- l d - U/2)/(2UT). 

With such an S-matrix we can readily determine the 
Bethe ansatz equations together with the solutions of 
these equations which populate the ground state. In the 
absence of a magnetic field, these are simply sets of real 
A and associated complex momenta, q±, given by 

g(Q a ±) = g(x(X a ) t iy(^a)) = K ± «/ 2 ; 

x(\) = U/2 + e d - VUT(X + v/A 2 + 1/4) 1 / 2 ; 



/ UT(-X+ V/A2 + 1/4) 1 / 2 . 



(2.58) 



This set of solutions is valid for + U/2 > 0. (But of 
course as e d = —e d — U, this corresponds to the condition, 
e d < —U/2, in which we are interested.) 

The above construction indicates that a hole scatter- 
ing at the Fermi surface in a system characterized by 
V,e d ,V p , and U is equivalent to a particle scattering at 
V = —V, e d = —e d — U,V P = —V p , and U. This is exactly 



what is to be expected from a particle-hole transforma- 
tion. And thus at the point e d = —U /2 = — e d , we obtain 
a second inequivalent description of scattering. Which 
description to use and when will be discussed in Section 
III. 



III. LINEAR RESPONSE CONDUCTANCE AT 

T = 

Here in this section we will compute the linear response 
conductance, G, at T = 0. To compute G we employ the 
Landaucr-Buttikcr formula, 



G 



■(T 2 (e = e F )+T?(e = e F )). 



(3.1) 



|0), (2.55) Here T is the transmission amplitude for electron scat- 
tering at the Fermi surface while T is related to ni mp , the 



number of electrons displaced by the impurity, through 
the relationship 



rri2 

i T/l 



2 At/i 



(3.2) 



(if Vl = Vn). Here 5 y^, the scattering of electrons off 
the impurity in the odd sector is given trivially by 5 ^ n = 
2 tan -1 (V p ). In the even sector on the other hand, n is 
given by <5 e j = 27rn imp | n — 2 tan -1 (V p ) as developed in 
the previous section. We once again emphasize that n imp 
is not necessarily equal to the number, n d , of electrons 
sitting on the dot but may include contributions from 
displaced electrons in the leads. Langreth makes this 
precise point in his proof of the Friedel sum rule^S As 
we will see this point will be crucial in interpreting our 
results. 

In the presence of finite V p , the structure of G is highly 
non-trivial. This is readily seen. Let us suppose that 
U/2 + e d > 0. In this case, the number of electrons 
displaced by the impurity is given by 



—D 



rQ fB 
nimpt = / G?A<7 imp (A) + I dqp imp (q); 
JQ 
rQ 

^im P l = / cL\<7 imp (A). 
JQ 



(3.3) 



The integral equations controlling <7i mp //9i mp are 



Pim P (g) = A(q)+g'(q) / dXa^X - g(q))a lmp (X); 
JQ 

rQ 

a imp (A) = A(A) — / dA'a 2 (A-A')cr imp (A) 
JQ 

+ / dqax(X- g(q))p imp (q). (3.4) 

J-D 

These equations cannot in general be solved in closed 
form. However examining the source terms, A.(q) and 
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(q- £ d ) 



FIG. 7: Plots of A(q) for a number of different values of V f 
Here T = .1 



at V p ~ 0.2 for the particular values of parameters cho- 
sen). At V p = V^ cnt + + , the S- function reappears with 
weight -1. Additional increases in V p serve to broaden 
and reduce the peak. 

To study further the linear response conductance (i.e. 
the behaviour of the equations (|3.4|) 'l we consider the 
cases of H = and H ^ separately. We first take 
up the zero field case. 

A. Zero Field Linear Response Conductance 

We will divide our discussion into three cases: 1) 
U + 2e d > 0, 2) U + 2e d < 0, and 3) U + 2e d = 0. 
The first two regions are studied through the two differ- 
ent developments of the system's ground state found in 
Section II (i.e. in the first 'particle' basis was 

used to develop the Bethe ansatz equations while in the 
second case we used a 'hole' basis). In the third case, 
U + 2e d = 0, the two developments are both ostensibly 
valid but give different results. To distinguish which so- 
lution should be employed, we examine the energy of the 
two corresponding ground states, choosing the state with 
the lower energy. 

1. U + 2e d >0: 



A(A), leads to significant insight. In Figurc[7]are found 
plots of A(q) for several values of V p . Recall A(q) is re- 
lated to the bare scattering phase of single particle eigen- 
states with energy q via A(q) = d q A(q)/2ir. A(q) is 
continuous in V p . As V p is made finite, the particle-hole 
breaking due to V p both shifts the peak away from q = e d 
(the direction of the shift is dependent upon the sign of 
V p ) as well as making the peak narrower. This is reflected 
in the general analytic form for A(q), 



- ' v{q-e d + V p T) + f 2 ' 

where f = 17(1 + V^ 2 ). 

While A(q) is continuous with variations in V p , A(A) 
is not. A(A) is given in terms of the bare scattering 
phase, not of single particle eigenstates, but of two- 
particle bound states. We plot A(A) in Figure |H1 for a 
range of values of V p . Two marked transitions are appar- 
ent in this plot. 

As V p is reduced from — oo a peak in A(A) develops. 
With V p approaching this peak develops into a 6- 
function of weight — 1. At V p = 0, the 5-function dis- 
appears and A(A) is smooth. As V p enters the region 
of positive values the delta function reappears but with 
weight 1. As V p is increased from + this <5-function at 
first broadens out (see V p = .1 in Figure |SJ but then 
narrows again. At some finite value of V p , designated 
Vp Crit , the 8— function disappears (in Figurc[H]this occurs 



a) U + 2e d = 0+ 

As a starting point, we begin by computing the con- 
ductance at U + 2e d = + , that is, at a value of the dot 
chemical potential just in excess of the point at which 
the model would possess particle-hole symmetry in the 
absence of a finite V p . At U + 2e d = + , Q = — oo, just 
as with the V p = Anderson model, and it is straight- 
forward to show that, 

/OO -j^ poo 

d\a imp {\) = - / dXA(X,V p ). (3.5) 
-oo * J — oo 

The discontinuities in A(A) as a function of V p just 
discussed are then intimately related to discontinuous 
changes in nj mp as a function of V p . We plot nj mp = 
"dmpT + n impi as function of V p in Figure OH We see for 
V p negative, rt; mp = 0. As V p is tuned to exactly zero, 
"imp = 1- Upon further increasing V p , nj mp changes to 
2. This marks the first transition in A(A) as V p is varied 
from — oo to +oo. When V p is tuned through its second 
transition, nj mp jumps back to 1 and then to 0. Given 
that we have 

e 2 

Gyi = ^sinV^impTA - 2tan~ 1 (Vp)), 
G experiences a similar set of jumps as shown in Figure 

nni 
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FIG. 8: Plots of FA(g) for a number of different values of V p . Here we have chosen U — 1 and T = .1. 



The jump in 7ii mp and G as V p transits V p = is the 
easiest to understand. As V p is changed away from 0, 
the particle-hole symmetry present is broken. With this 
symmetry present, nj mp necessarily is 1. This change 
by ±1 in rii mp is mimicked by an exact diagonalization 
calculation of the ground state occupancy for a lattice 
model equivalent to the continuum Anderson model. To 
be specific, we consider the lattice model 



H 



lattii 



■* (4a C i~^ + h.C.) 



i=-^- + l,a 
(JV-l)/2-l 

-* (cLci + i CT + h. c .) 

+ V ^2( c L c o<r + cl lo .c 0<7 + h.c.) 

(7 

+ V LR J2(cLc~la+h.C.) 
a 

+£d^2cl a c 0a + C/c^cotc^cox- (3.6) 

On the lattice, the dot sits at i = 0. Here t is the hopping 
amplitude between the "lead" (i.e. non-dot) lattice sites 
and is roughly equivalent to the bandwidth. V is the 
hopping amplitude between the dot site and the neigh- 



boring lead sites, i = ±1. The hopping amplitude, Vlr, 
allows the electron to bypass the dot, hopping from i = 1 
to i = — 1 and vice versa. The dot parameters U and €d 
are the same as in the continuum Hamiltonian 12.1fl . 

We have constructed the ground state of this lattice 
model for a number of small lattice sites N < 13, TV odd, 
at the point U + 2ed = 0. With Vlr = the ground state 
is degenerate, one state with N — 1 electrons, one state 
with TV + 1 electrons. Averaging over these states, each 
site on the lattice contains exactly one electron. Turning 
on Vlr breaks the degeneracy of these two states. For 
Vlr > 0, the system prefers to be in a state with TV — 1 
electrons, while for Vlr < the system prefers to be in a 
state with TV+ 1 electrons. Thus as Vlr is swept through 
the ground state occupancy changes by 2. We suggest 
this is the lattice equivalent of the change by 2 as seen in 
n; mp as computed by the Bethc ansatz. 

Now we point out, as is obvious from Fig. [5] and 
Fig. ^2 that the changes in nj mp and n™ as V v passes 
through zero are opposite in sign. In the context of the 
Bethe ansatz computation this is readily understandable. 
Consider V p negative. The exact diagonalization com- 
putation indicates for Vlr negative the system should 
increase its occupancy by one electron. For the Bethe 
ansatz to increase the occupancy of the system, it must 
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FIG. 9: A plot of the total number of electrons (both spin 
species) displaced by the presence of the impurity as a func- 
tion of V p . We have assumed U + 2ed = + and have chosen 
U — 1 and F = .1. The plotted filled circles indicate that at 



V v = and V v = K c 



0.19, nimp is precisely 1. 



FIG. 10: A plot of the conductance through the impurity as a 
function of Vp. We again take U + 2ed — + as well as U = 1 
and r = .1. The plotted filled circles indicate the values of G 
V p crit ~ 0.19. 



add a two-particle bound state to the bulk distribution. 
For an overall change of one electron, the impurity oc- 
cupancy must then compensate by decreasing by one. A 
similar argument applies to V p positive. The opposing 
changes in ni mp and n™ arc then ultimately a conse- 
quence of the Bcthe ansatz constructing the ground state 
out of two-particle bound states. 

This change in lattice occupancy can be characterized 
as a many body effect. If N = 2R + 1 with R even, 
the change in ground state occupancy is present only for 
finite U. If U = instead, the total ground state occu- 
pancy is precisely N regardless of the value of Vlr (see 
Figure 11). R must be even for only then does the non- 
interacting lattice system have zero energy single particle 
excitations. With R odd (and no zero energy single parti- 
cle excitation), the ground state occupancy jumps as Vlr 
passes through zero regardless of the value of U. But in 
choosing R for taking the continuum limit, R even is ap- 
propriate as the continuum limit possesses such a zero 
energy excitation. 

Although the changes in the ground state occupancy 
were computed for small lattice sizes, the nature of the 
results are robust and will persist to arbitrarily large lat- 
tices. 

The origin of the second transition in n,j mp occurring 
at finite positive V p is less transparent. The transition 
at V p = 0, on the basis of the exact diagonalization cal- 
culations, can be interpreted as one electron on average 
being added or removed from the ground state. However 
at the second transition no obvious change in the ground 



state occupancy occurs in the exact diagonalization of 
the lattice. A clue to the nature of the second transition 
can be found in the value of V p at which it occurs. 

To determine this value of V p we need to first determine 
the value of A at which A (A) develops a singularity when 
V p passes through the second transition. At this value of 
A, the argument of the inverse tangent appearing in the 
scattering phase of the two-particle bound state, i.e. 



^bd. state (A) = 

satisfies 

Im 



-2Rctan~ 



r 



_x(A) - iy(X) - e d 



v„ 



= 1. 



a; (A) - iy(X) - e d 
Solving this equation for large U/T we find 



(3.7) 



A 



u r 

8T 2U' 



The V„ = singularity appears at A = ^ — ^ while 

To 



the V p = Vp rlt singularity is found at A 



sr 2U 
sr ^ 2U- 



determine then V^ cnt , we insist that the argument of the 
inverse tangent in <5bd. state satisfies 



Re 



F 



Vp 



_x(X) - iy(X) - e d 
Solving we find the critical value of V p is 



= 0. 



V, 



crit 



2r 
77 
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b) U/2 + Ed finite and positive: 
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FIG. 11: Plots as a function of Vlr of the deviation in the 
ground state occupancy of an N-site Anderson lattice model 
from a state with occupancy N—l. The exact diagonalization 
calculation of the model takes U + 2ea — 0. The plotted filled 



circle indicates that at Vlr = 0, An a 



1 corresponding to 



exactly one electron living on the dot site i = 0. 



We now move on to consider how the linear response 
conductance varies as Ed is increased from the value 
—U/2 + + . G evolves according to the value of nj mp f n 
at Ed = — C//2 + 0+. If n- lmp y^ = (as will be the case for 
V p < as well as for V p > V^ crlt ), as Ed tends to large val- 
ues, the conductance, G, will return to the value seen at 
Ed = —U/2 + + . We have that G depends upon nj mp j n 



G = 2— sin 2 (7m impTA - 2tan- 1 (V p )). (3.9) 



nmp 



j/l is necessarily ~ at Ed/T 1 for in this limit 
the Fermi surface, Q, of the A-excitations tends to oo and 
so 



lim 



"impT/l 



= lim 



dACTimp(A) = 0. 



Thus if riimpT/l 
have 



at Ed = —U/2 + + , we necessarily 



G(E d = -17/2 + 0+) = G(E d = oo) 



Thus the value of V p at which the second transition oc- 
curs is roughly 2T/U. To interpret this we suppose that 
at large U/T, the Coulomb repulsion U serves as a rough 
cutoff for the theory. With V p a dimensionless parame- 
ter, the energy scale corresponding to V^ crlt is then 2r or 
roughly the level broadening of the dot. Thus the second 
transition occurs when the tunneling from one lead to the 
other induced by V p has grown to be the same magnitude 
as tunneling through the dot. 

This transition only occurs at V p positive. At V P U ~ T, 
the potential scattering term is strong enough to remove 
the electron from the impurity itself. That is, the elec- 
tron that must be ejected from the system (at V p > as 
indicated by the exact diagonalization computation) is 
now removed directly from impurity. We no longer have 
the construction by which V p > induces the removal of 
a two-electron bound state in the bulk and a concomitant 
addition of an electron to the integrated impurity density 
of states. 

For U/2 + Ed > 0, this suggests an inequivalcnt po- 
sition for particle vs. hole scattering. There is no cor- 
responding transition at V P U ~ — T where the impurity 
occupancy is increased by one. We however find such 
a transition for U/2 + Ed < where the roles of parti- 
cles and holes are reversed. Exactly at U/2 + Ed = 0, 
both transitions, V P U ~ ±r appear, as they must by the 
demands of a particle-hole transformation which leaves 
U/2 + Ed invariant but changes the sign of V p . 



2— sin 2 (2 tan" 1 (Vp)) 
h 



4Kf 



h (1 + Vp 2 ) 2 ' 



(3.10) 



The value of G at these two values of Ed is then solely a 
function of V p . 

There will of course be significant variation in G as 
Ed changes from —U/2 + + to +oo. This variation will 
occur at the value of Ed at which the A-Fermi surface, 
Q, reaches the position of the largest variation in A(A). 
As V p is tuned away from 0, the position at which the 
non-analyticity in A (A) opens up is A ~ U /8T. It is thus 
reasonable to locate the largest variation in A(A) at this 
value of A. The behaviour of Q as a function of Ed has 
been determined in Ref. |3(J Provided Q > 0, Q can be 
written in terms of Ed, U , and T as follows 



Q 



<f 



l 

2^ 

U/2 



log(2?reg*); 



V2T7T 



(3.11) 



Thus the critical value of Ed at which Q = U/8T and G 
sees a large variation is 



U 
~2 



U U 

r 25/2 



provided U/T > 1. 
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The size of the variation in G occurring at e d nt can be 
estimated. From (|3.9[1 . the variation in G can be related 
to the variation in 



limpj 



limp J, • 



SG 



9 eV 9 Vpjl-V?) : 

2\2 V u "imp^ 



(1 + ^p 2 ) 



(3.13) 



i5ni mp can in turn be estimated to be 



M Re( 



V„ 



7T i + y p 2 
i i 

7T 1 + F p 2 ' 



'x(A) - iy(A) - e<i' 



(3.14) 



The first estimate in (|3.14|l follows as we expect 
obeys 



''imp 



/•Q /-Q 
iim P = 2 / rfA(7 imp (A) = / d\A(X) 
Jq Jq 

= -Rc5{x(X)+iy(X)). (3.15) 



The second equality in l|3.14|l follows provided V p is large 
while the third and final estimate holds as Re(r/(x(A) — 
iy(X) — €d)) varies between —1/2 and 1/2 in the vicinity 
of A = Q. Hence SG is given by 



SG 



^p(i-K 2 ) 



1 - 3V P 2 



(1 



K 2 ) 3 



(i + K 2 ) 4 



(3.16) 



According to this expression, near V p = 0, the sign of the 
variation changes. 

In Figure 12 these various analytic formulae are borne 
out. At V p = —0.5,-0.1,0.3, and 0.75, we have 
n imp |/|(C//2 + td = + ) = 0. We see that the value 
G asymptotes to as — > oo equals the value of G at 
U/2 + Ed = + . The variation in G takes place at 
(_d ~ —0.1 in accordance with l|3.12|l . Finally the size 
of the variation in G clearly depends upon V p with the 
variation decreasing in size as V p increases. We also see 
that variations in G on opposite sides of V p = are of 
different signs in accordance with (|3.1()|) . 

We now consider the behaviour of G as a function of 
e d in the region of V p , < V p < V^ Tlt , where n impT/i (e d + 
U/2 = 0+) = 1. By virtue of we again find that 

G(e d = -U/2 + Q+) = G(e d = oo) 
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FIG. 13: A plot of the total number of holes (of both spin 
species) displaced by the presence of the impurity as a func- 



FIG. 14: A plot of the conductance through the impurity as 
a function of V p . Here U + 2ed = CP while U — 1 and F — .1. 



tion of V p . We have taken U + 2td = , U = 1 and The plotted filled circles indicate the values of G at V p — 



.1. The plotted filled circles indicate that at V p 



and V p = Vp= Tlt 



-0.19, n^p is precisely 1. 







h (l + V 2 ) 2 ' 



(3.17) 



The variation in G as e<j is increased from —U/2 still 
takes place at = e d Tlt (|3.12|) as the primary variation 
in A(A) still occurs at A = U/8T. However the size of the 
variation in G differs from the case ni mp yi(U/2 + ed = 
+ ) = 0. Here G will take on its maximal (2e 2 /h) and 
minimal (0e 2 //i) values as ed is increased upwards from 
—U/2. This is a consequence of <|3.9ll and the fact that 
n impT /I wm " decrease monotonically from 1 to as e d is 
made large and positive. 

This behaviour in again seen in the conductance plot 
for V p = 0.1 in Figure 12. G there is 0.08e 2 //i at e d = 
—U/2 + + , varies rapidly in the region of ed = e" 1 ' = 
—0.1, taking on values of both 2e 2 /h and 0e 2 /h, then 
proceeds to decay back to 0.08e 2 /h for larger values of 
ed- 
Having considered values of V p where nj mp |/|(e f / = 
— U/2 + Q + ) = and 1, we now turn to the two remaining 
undiscussed values of V p , V p = and V^ crlt , where here 
n impT/l( £ d = —U/2 + + ) = 1/2. For these two values of 
V p , the two limiting values of G are 

2e 2 



G{e d = -U/2 + Q+) 



G(e d = oo) 



h 



cos 2 (2tan- 1 (y p )) 



= 2 



1-K 



h V 1 

2e 2 



V 2 

AV 2 

p 



h (1 + V 2 ) 2 ' 



(3.18) 



For V p = and V^ crlt , A(A) does not possess large varia- 
tions at A = U/8T. Thus G sees no large variations. At 



Vn 



0, G decreases monotonically from 2e jh as td is 



increased from —U/2. If U/T 3> 1 and correspondingly 
^/crit ^ 2 5 G(V P — Vp Tlt ) will also decrease monotoni- 
cally as ed is varied upwards from —U/2. The plots of G 
for V p = and Vp Crlt in Figure 12 bear this out. 



2. U + 2e d < 

To treat the case U + 2td < 0, we need to use the 
development of the Bethe ansatz equations detailed in 
Section II. F. Here everything is cast in terms of holes 
not particles. The conductance is given by 

G(e d , V p ) = 2^ sin 2 ( 7 rn| 1 n 1 1 p c T/i (e d , V p ) + 2t^- 1 (V p )) 

(3.19) 

where ^ifit/i is the number of holes displaced by the 



'impt/i 
.j hole 
t impT/l 



impurity. ^j^DT/i ^ s gi yen by (at zero field) by 



n hole 



Q 



n impT/i = / dXcrimpiX, -V p ,-U- e d ), (3.20) 
Jq 

where CTi mp (A, — V p , — U — ed) is governed by 13.4fl . 
a) U/2 + e d = 



Akin to the treatment of U +2ed > in Section III.A.l, 
we first examine the case U + 2ed — _ . At U + 2e d = _ , 
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the number of holes displaced by the impurity is 

1 r°° 

<; m (U + 2e d = 0- V p ) = - d\A(X 1 -V p ). 

L J -oo 

(3.21) 

Thus n™; yi (U + 2e d = 0",F P ) = n impUl (U + 2e d = 
+ , — Vp) as would be expected from a particle- hole trans- 
formation. Plotted in Figure 13 is n^(U + 2e d = 0") 
as a function of V p . It is a mirror image of Figure 9. 

We ascribe the jumps in nf°^(U + 2e d = 0~,V P ) at 
V p = and V p = —V p mt to similar mechanisms as 
the jumps seen in n,i mp (U + 2e d = + ). The jump in 
n tmp(U + 2td = _ ) as V p moves away from is caused 
by the breaking of particle-hole symmetry by finite V p . 
For Vp = , the system wants to have one more/one 
less hole in the ground state. It thus adds/subtracts one 
two- hole bound state to/from the bulk occupancy and 
subtracts/adds one hole from/to the impurity occupancy. 

Again the jump in n\^° at V p = — Vp Crlt can be thought 
to be due to the rough equivalence occurring at this point 
of the two tunneling paths through the impurity. Here 
the transition occurs at V p negative. The hole that must 
be removed from the system at V p < now is removed 
directly from the impurity occupancy as opposed to being 



removed from the bulk. 

As in the case U/2 + e d = 0~, the conductance at 
U/2 + e d = 0~ experiences a set of discontinuities. This 
is plotted in Figure 14, a mirror image of Figure 10. 

b) U/2 + e d finite and negative: 

The behaviour of G as e d is decreased from —U/2 is 
determined by the value of n^^y^ at e d — —U/2 + CP. 
As we have already carried out an analogous discussion 
of G's evolution as e d is increased from —U/2 and its 
dependence upon rii mp y^(e d = —U/2 + + ), we merely 
summarize in point form the results in this case. 

a ) <p T /l( e d = -U/2 + 0-) = 0, (V p > and V p < 
_ V crit ): 

• The conductance at e d = —U/2 + CP equals G at e d — 
—oo, namely 

AV 2 

G(e d - U/2 + 0-) = G(e d = -oo) = (i - ^ 2)2 . (3.22) 
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• The largest variation in G as ed varies between — oo to 
-U/2 + (T will occur at 



-rcrit 



-U-ef* 



u u u 
1 ~ V r 2572 



^ l0g(27re 8f ) 



The variation in G at e<j = e™' is given by 

v p (i-v 2 ) i- 3 y p 2 + ^ 4 



(5G = 2 — 
h 



(i + K 2 ) 2 



(i + ^ 2 ) 4 



(3.23) 



(3.24) 



• These features are seen in Figure 15 in the plots of G 
for V p = .75, .3, .19, .1, and -.5 (where lf lt ~ -0.9). 
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b ) n | 1 mpT/l( £ d = -U/2 + 0-) = 0, (0 > V p > -V p rit ): 

• At the two limits, = —U/2 + CP and = -co, the 
conductance satisfies 



G{e d - U/2 + 0") = G(e d = -oo) = 2 



e 2 

h {i+v^r 

(3.25) 



crit 



• The variation in G occurs again at e" = —U — e d ] 

• The variation in G as ed is varied from — oo to — U/2 + 
0~ is 2e 2 /h. That is, the minimal value of G is 0e 2 /h 
while its maximal value is 2e 2 /h. 



FIG. 16: A plot of the impurity energy at U + 2e<i = as 
a function of V p with (7=1 and T = 0.1 as computed with 
both the particle and hole Bethe ansatz solutions. 



s 2 1 



T^crit q 



6Vp Crit " 



[1 + 



(3.27) 



A similar zero temperature discontinuity in G is seen at 

stion 
17/2. 



V p = — Vp rlt . This discontinuity gives rise to the question 
of what is the actual value of G when ed is precisely 
We resolve this question in the next section. 



• These features are seen in Figure 15 for the conductance 
plot for V p = -0.1. 



U + 2e d = 



c ) <p T /l( e d = -U/2 + 0-) = 1/2, (V p = 0, -V^): 

• At the two values of the dot chemical potential, td 
—U/2 + 0~ and = — oo, the conductance satisfies 



G(e d -U/2 + Q-) = 2 



e 2 (1-V P 2 ) 
h 1 + V 2 ' 



G(e d = -oo) = 2- 



h (1 + V 2 f 



(3.26) 



• G varies smoothly and monotonically as is varied 
between — oo and —U /2 + 0~. 

• Comparing Figure 12 and Figure 15, we see the con- 
ductance is continuous as td crosses —U/2 for all values 
of V p but for Vp = Vp 11 *. This jump is of magnitude 

AG = G{e d = -U/2 + 0+) - G(e d = -U/2 + 0~) 



At e d = —U/2 (and td = —U/2 alone), two different 
developments of the Bethe ansatz equations are avail- 
able: one involving a basis of single particle states and 
one involving a basis of single hole states. These two 
treatments are not equivalent for all values of V p as the 
discontinuity at e d = —U/2 in G at V p = iV^"* shows. 
This discontinuity is a reflection of our ability to con- 
struct two not one putative ground states. To determine 
the true ground state at a given V p , we compute the en- 
ergy of the two competing states. 

The energy of the particle Bethe ansatz state (for H = 
0) equals 

f Q 

E palt (e d ,V p ) = L dXa(X,V p ,e d )2x(X), (3.28) 

JQ 

where cr(A), the density of bound states at a given A, 
satisfies (|2.35f) and 2a;(A) is the (bare) energy of a two- 
electron bound state. Decomposing <r(A) into its bulk 
and impurity contributions, a = Ubuik + &vmp/L, allows 
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us to similarly decompose -Epai-t^ 



£part(ed, V p ) - Sp" r t (erf, Vp) + —-Bp 



imp _ 

art 5 



Kar P t( e d,^>) = / dXa imp {X,V p ,e d )2x(X,e d ,U); 



E^(e d ,V p ) = L J dXa halk (X,e d )2x(X,e d ,U). 

(3.29) 

(Tb u ik(A) shares its independence from V p with E^%£ . The 
energy of the hole Bethe ansatz state is given by 



S h oie(ed, Vp) = L J dXa holc (X, V p ,e d )2x(X, -U - e d ,U) 
+V P + U + 2e d 

= L j dXa(X,-V p ,-U-e d ) 
x2x(X,-U-e d ,U) 

+V P + U + 2e d . (3.30) 

We have used the fact that the density of two- hole bound 
states in the hole Bethe ansatz state, <7hoie(A, Vp, e d ), is 
equal to the density. cr(A, — V p , —U — e d ), of two-electron 
bound states as computed using the particle Bethe ansatz 
up to a change in sign of — V p . We can again perform a 
decomposition into bulk and impurity pieces: 

Eh \ c (e d} Vp) = E^(e d ,V p ) + —E 1 ™?; 

*%Z{e d ,V p ) = V p + U + 2e d 

+ [ dXa imp {X,-U - e d ,-V p )2x(X,-U - e d ,U) 



E hole ( e d, Vp) 



L / dXa bulk (\ -U - e d )2x(X, -U - e d , U). (3.31) 




FIG. 17: A plot of the conductance as a function of V p for 
U + 2e d = 0, U = 1 and T = 0.1. The filled circles mark 



ycnt 



±0.19. 



symmetric, i.e. 

G{e d = -U/2,V P ) - G{e d = -U/2,-V p ), 
as it must be. 

The determination of the true ground state as a func- 
tion of Vp thus eliminates a potential problem. If one 
were to look at G(e d = —U/2) on the basis of either the 
particle Bethe ansatz state or the hole Bethe ansatz state, 
G would not satisfy particle-hole symmetry (see for ex- 
ample Figures 10 and 14). Only by taking into account 
both solutions, do we establish a conductance consistent 
with particle-hole symmetry. 



So while and E$J% are equal at e d = -U/2, the 

impurity contributions are not. By determining these 
impurity contributions then we can determine the state 
of lowest energy, the true ground state. 

Plotted in Figure 16 are the impurity contributions to 
the energy arising from both the hole and particle Bethe 



imp 

part, hole 



IS 



ansatz states, c^mp, necessary to compute E 
determined by equations (|2.37|l and Q2.41fl . For V v > 0, 
we see the particle Bethe ansatz state forms the true 
ground state while for V p < 0, it is the hole Bethe ansatz 
state that gives the ground state. At V p = 0, the two 
states become equivalent. 

Having determined the true ground state as a func- 
tion of V p , we can determine G at U = 2e d . We plot 
G(e d — —U/2) in Figure 17. As a function of V p , it 
sees three discontinuities: two at V p = ±V P and one at 
V p = 0. Despite these discontinuities, G is particle-hole 



B. Finite Field Conductance 

In this section we will consider how G behaves as a 
function of H . Making H finite changes the ground state 
of the system. From consisting solely of two-electron (or 
hole) bound states, the ground state begins to admit the 
presence of single particle states with real wavevectors, q. 
The scattering phase is still given by (|2.47|) and (|2.51|) but 
now the number of electrons displaced by the impurity is 
given by (|23U|l with B > —D: 



"impT = / dAcr imp (A) + / dqp imp (q); 
iQ J-D 



"impi = / dXa imp (X). 
iQ 



(3.32) 
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FIG. 18: Plots ol the impurity density ol states pim P /d q e as a lunction of energy for a number of a values of V p . Here we have 
taken [7 = 1 and T = 0.1. 



Thus the behaviour of Pimp(<z) is dispositive of how G 
will vary as H changes from zero (as H is increased, B 
increases from —D). 

To obtain a feel for how G evolves as a function of H, 
we examine Pi mp (q) at both H = and Ed = — U/2 + + . 
If there is low energy spectral weight in pi mp (q), we ex- 
pect G to vary rapidly with H. If not, G will be relatively 
constant as a function of H. We plot Pim P (q) / d q e(q) vs 
e(q) for a number of representative values of V p in Figure 
18. e(q) is the energy needed to add a g-excitation to the 
system. As there are no g-excitations in the groundstate 
at H = 0, e(q) is necessarily positive. The derivation of 
e(q) may be found in Section II of Ref. |2f| (applicable to 
the case here as e(q) does not depend on V p ). 

The behaviour of Pim P (<?) at H = and U/2 + Ed = + 
as a function of V p can be put into three categories. If 
n- wvVi {V p ,e d = -17/2 + 0+) = 1/2 (i.e. V p = 0, V^) 
there is low lying spectral weight. With V p — 0, this 
spectral weight is associated with (but not identical to) 
the Kondo peak in the dot density of states. This low- 
lying weight ensures G varies rapidly as a function of 
H. If n impyi (V p ,e d = -U/2 + 0+) = 0, i.e. V p < or 
V p > Vp Tlt , the low- lying spectral weight doubles approxi- 
mately in size. It is as if spectral weight has shifted from 
being below the Fermi energy and encoded in n imp -f 
at V p — to being above the Fermi energy and en- 
coded in pi mp (q). Again this spectral weight guaran- 
tees G will change rapidly with the introduction of fi- 
nite H. Finally if n imp1/l (V p ,e d = -U/2 + 0+) = 1, 
(i.e. < V p < V^ crlt ), the low- lying spectral weight in 
Pim P (q) is absent. Here it is as if the spectral weight 
has shifted from p- lmp (q) to below the Fermi surface in- 
creasing n imp yi(Vp, ed = —U/2 + + ) from 1/2 to 1. In 



this case G will be unresponsive to small magnetic fields 
(H ~ T K ). 

1. Conductance at td — — U/2 + + 

We first turn to computing the conductance at the 
point Ed = —U/2 + + . In Figure 19, we plot the con- 
ductance as a function of H for Ed = —U/2 + + . We see 
that there is significant variation in the conductance for 
V p outside the range (0,V^ crlt ). In contrast for V p inside 
this range, G(Ed = —U/2 + + ) is nearly constant. This 
is in accordance with the presence or absence of low-lying 
excitations as parameterized by Pi mp (q)- 

Figure 19 also encodes information on the finite field 
conductance for Ed = —U/2 + Taking V p — > — V p 
in Figure 19 gives the behaviour of G(E d — U/2 + 0~). 
We thus realize G experiences generically a discontinuous 
jump at finite V p and finite H as Ed crosses —U/2. We 
plot the size of this jump as a function of H for a set 
of representative values of V p > in Figure 20. The 
corresponding discontinuities for V p < can be obtained 
via 

AG(V P ) = AG(-Vp). 

While only at V p = V^ rit is AG non-zero at H = 0, 
as H is turned on, AG rapidly becomes finite and large. 
Although not plotted, these discontinuities will persist to 
large values of H and will not disappear until H yC/T. 

We are able to derive analytically some results at 
E d = -J7/2 + ± . We will focus upon E d = -U/2 + 0+ 
as results from E d = —U/2 + 0~ can be obtained via 
a particle-hole transformation, i.e. V p — > —V P . With 
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£d = — U/2 + + , A(A), the bare impurity density of 
states can be approximated via 

A(A, V p ) = A(A, V p = 0) + S(V P )S{X - A ), (3.33) 

provided |Vp| is not much greater than V^ crlt . If in ad- 
dition U r, this approximation works for the entire 
range, < V p < V" lt , as V^ crlt is small. We will work 
under this assumption. Then from our analysis of the 
H = conductance we can write 



A ~ U F (1 2 Vp )■ 
8r 2U y Vp= Hth 



o<v p < y; 1 



(3.34) 



and 



0, 

S(V P ) = { i. 



r -1, for V p = 0- 
for V p = 0; 



for < V p < Vp crit ; 



0, for V p 



-1, for V p = Vp crit 



(3.35) 



H 



In the expression for Ao, we have used a linear extrapo- 
lation in Vp between the value of Ao at V p — and Ao at 

V„ = V crit . 
* p v p 

To compute the magnetoconductance at H ^ and 
€d = —if/2, we need to solve ultimately the integral 
equations 



Pimp(g) 

Cimp(A) 



As the limits of the integrals involving <7i mp (A) run from 
— oo to +oo, the above equations can be recast into a 
single equation: 

/oo 
dXA(X)s(g(q) - A) 
-OO 

-g'(q) I dqR(g(q) - g(q'))p imp (q), (3.37) 



A(g)+fl'(g) / dX ai {g(q) - X)a imp (X); 

J — OO 

A(A)- / dAa 2 (A'-A) ( 7 im p(A') 

OO 

dqa x {\ - g(q))pim P (q)- (3.36) 



where 



S (A) 



2 cosh(7rA) ' 



R(X) = ± / cL,e^-i 



M 



(3.38) 



Needed information on <7i mp (A) can extracted from the 
relation 



1 f°° ~ If" 
dXa imp (X) = - / A(A) - - / Pimp (q). (3.39) 

i 1 J-oo z J-D 



To solve the integral equation (|3.37|) , we note that using 
(13.331) and the results borrowed from Ref. l36t the source 
term of i|3.37[) can be recast as 



A(<?) + g'(q) / d\A(X)s(g( q ) - A) 

J — OO 

= g'(q)(S(V p )s(g(q) - A ) - s(g(q) - J- 1 )). (3.40) 
where 

8r 2C 

and again we assume V p to be small. With this 
observation, the integral equation can be solved and 

J-d dqpim-p(.q) determined from Wiener-Hopf techniques 
on exactly the same lines as in Ref. |3(J The result is 



D 



dqPimp(q) 



V2 ~ G+(i7r(2n + l)) ^_ wa ^ 



7T <— ' (2ra+l) 

n=0 V ' 



x (-1)" ^ e -(2«+D//) _ S{V p )e^ 2n+1 ^j ; 

(3.41) 



■ i w 2 /err 

6 = -log — W ). 

7r H \ ne 



We point out that for £//T > 1, I^ 1 = A . Thus when 
< y p < V p crit and S(V P ) = 1,J B D dqp imp (q) ~ 0. 

For small H/Tk, where the Kondo temperature Tfc (at 
e d = -U/2) is defined by^S^ 



T k = 



(3.42) 



we can arrive at explicit expression for G(V p ,ed = 
-U/2 + 0+). With small H/T k , the sum in lpOT|l can 
be approximated by its first term. So using l|2.9)l . I|2.47|l . 
(|3.39|l , and (|3.41|) , we can write down the following 



G(V p ,e d = - 



+ ,^«l) 

J-k 



e 2 (tt 2 H 2 l + K_-6V£ 

L h \~m (i+v^r 

16 T k J ' 



'1 



4K 



for v p < o, v p > y p crit + o+ 

for Vp = 0; 

for < Vp < y p crit ; 

for V p = I/ p crit . 



(3.43) 
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FIG. 19: Plots ol the magneto-conductance at e d = —U/2 as a function of H for a number of a values of V p . In these plots, we 
have taken [7 = 1 and V = 0.1. 
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FIG. 20: Plots of the discontinuity in the magnetoconduc- 
tance at e d = — U/2 as a function of H for a number of a 
values of V p . Here the discontinuity is defined as AG(H) — 
G(e d = -U/2 + 0+,H) - G{e d = -U/2 + 0',H). In these 
plots, we have taken [7 = 1 and F — 0.1. 

We see the deviations in G from its zero field value are 
quadratic in H/Tk. This indicates that regardless of the 
value of V p , the ground state is Fermi liquid-like at = 
—U/2 + + . Thus the discontinuous changes as V p is 
varied mark first order transitions between Fermi liquid 



states. 

On the basis of these expressions, we are able to write 
down analytic expressions for the discontinuity in the 
conductance occurring at small but finite H and V p 
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FIG. 21: Plots ol the V p = 0,0.1,0.19, and 1 magnetoconductance as a function of td for a number of values of H. In these 
plots, we have taken U = 1 and F = 0.1. 
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AG = G(V p ,e d = ~+0+)-G(V p ,e d = 



G(V p ,e d = ~ 



4+»-> 

.Q+)-G(-V p ,e d = ~+0+) 



— 

u h ' 



8V„ 2 



J_ 

16 



ft, ' 



4 (i+v p 2 ) 2 



for Vp = 0; 

for < V p < y p crit ; 

for V v = K crit ; 



for y p = V p 



crit 



(3.44) 



For < V p < V p CTlt , the jump in the conductance devel- 
ops as H 2 /T 2 from 0. At V p = V p crit the magnetic field 
decreases the size of the jump. At Vp = V p rlt + + the 
discontinuity disappears for small H/Tk- But for larger 
values of V p , as evidenced by Figure 20, the jump reap- 
pears. 



2. Finite H Linear Response Conductance Away from 
e d = -U/2 

Away from e d = —U/2, the transmission amplitude 
of the electrons is determined by (|2.47|) where now Q ^ 
— oo. In general these equations do not admit analytic 
solutions. Nonetheless we can solve them numerically. In 
Figure 21 we plot the linear response conductance for a 
variety of values of H and V p as a function of e d . We take 
V p > (again V p < can be obtained via a particle-hole 
transformation, G(e d , V p ) = G( — U/2 — e d , —V p )). 

In Figure 20a as a reference we have plotted G as a 
function of e d at V p = for successive values of H. We 
see that each plot is symmetric about e d = —U/2 (as it 
must be for V p = by particle-hole symmetry). As H 
is increased, the conductance peak at e d = —U/2 splits 
into two smaller peaks. For large H, these peaks occur 
at e d = H/2 and e d = —U — H/2 with a conductance 
G ~ e 2 /h. The peaks may be modeled by Breit-Wigncr 
resonances, i.e. 



G(e d ,H->T k ) = Lf sin 2( tan -i ( _[_ )) 



+ sin 2 ( tan 1 ( 



e d + U + H/2 



)) 



Moving to finite V p , the behaviour of G as a function 
of e d radically changes. No longer symmetric about e d = 
—U/2, G now sees sharp variations. These variations 
arise for the same reason variations at H — happen. 
They thus occur at gate voltages corresponding to 



€ d = £ d 
and 



-U/2- 



U U 

F 25/2 



^ l0g(27re 8f } 



_ U cr it 
~ e d ■ 



(3.46) 



These variations are absent only for V p — V p crlt , e d > 



-U/2 and V p = -V^' l \e d < -U/2. As discussed pre- 
viously, V p = Vp rlt is the only non-zero value of V p at 
which A(A, V p ) is relatively smoothly varying. 

From Figure 21, we see the discontinuities in the mag- 
netoconductance at e d = — U/2 discussed previously. As 
H is increased, the discontinuities grow smaller. In fact 
one can check (numerically at least) the discontinuities 
disappear in the regime H 3> V. This suggests that 
large H destroys all non-perturbative physics. Indeed 
for H > T, G(e d , H > T k ) behaves as 

G(e d ,ff»T fc ) = 
& - ( sin 2 ( tan"^— i— + V p ) + tan" 1 (V p )) 

r 



+ sin 2 ( t^- 1 ( ed + x /2 + u + V p ) + tan" 1 ^))) , 



(3.47) 



in analogy with 



h (e d - H/2) 2 + T 2 

r 2 

h (e d + U + H/2) 2 + T 2 ' 



(3.45) 



Thus G appears as if due to two species of spinlcss 
fermions, one with a resonant level at e d = H/2 and 
one with a resonant level at e d = —H/2 — U. 



IV. ATTRACTIVE U ANALYSIS 

In this section we will consider the U < Anderson 
model. In doing so we will demonstrate that the non- 
analyticities present at U > and V v ^ disappear. 
This indicates the problem is non-perturbative in U . 
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The Bethe ansatz equations describing the quantiza- 
tion of momentum are the same as in <|2.33[) . However 
for U < 0, the solutions of the equations are of a different 
character. At U > 0, the ground state consists of bound 
states characterized by complex pairs of wavevectors, q±, 
tied to particular values of A (the spin momenta). But 
at U < 0, the ground state is composed of real values of 
q independent (nominally) of any specific A. As in the 
U > case, we can write down a set of equations describ- 
ing the distributions, p(q) and <r(A), of fc's and A's in the 
continuum limit of the ground state: 



p(q) 



a(X) 



1 

2^ 



A(«) 



L 

Q 



g'(q) / dAai(A - g(q))a(X); 



dX'a(X')a 2 (X-X') 



dqp(q)a 1 (X-g(q)). (4.1) 



-D 



As with U > 0, the limits — D and Q mark the band- 
widths of the k and the A-distributions. But now Q oc- 
curs at A << - in fact it can be shown to be given by 
solving 



\u\/2 + e d - VW(-Q + (Q 2 + 1/4) 1/2 ) 1/2 



D. 



The remaining two limits, B and Q, mark the Fermi sur- 
faces of the two distributions and so are characterized 
by 

r-B 

dqp(q) 



D + p 



-D 

Q 



d\a(X) 



Q 



TT 

D + p 
2tt 



S z 



D + p H 



2- 



2tt 



(4.2) 



Here p is the chemical potential. We introduce p for 
technical reasons. If one wants to explore the effects of 
varying e^, one can alternatively vary p keeping fixed 
as it is only the relative distance between the dot energy 
level and the lead-electron Fermi surface that is relevant 
to the physics. As such for the purposes of any actual 
computations, we fix at — U/2 and simply vary p. 

As before we can divide these densities into a bulk and 
an impurity piece via writing p(q) = /Obuik(s) +Pimp(<z)/£ 
and er(A) = er bu ik(A) + o- imp (A)/L: 

1 f Q 

Pbuik(g) = 7T+9{(l) I d\ai(X~ g(q))a bu i k (X); 

27r Jq 



O bulk (A) 



dA'cr bu i k (A')a2(A - A') 



dqpbu\k(q)ai(X - g(q)), 



(4.3) 



-D 



and 



Q 



Pimp(g) = A(q)+g'(q) / dX(n{X - g(q))a imp (X); 
JQ 



o-imp(A) = - dX'ai mp (X')a 2 (X - A') 

B 

+ / dqp imp (q)ax{X- g{q)). (4.4) 

J-D 

The impurity spin and occupancy are given by 



«im P T + ™impi = / dqp imp (q); 
J-D 

^-impt T-impJ. = -^imp 

. rB 

= -2 / dXa imp (X) + / dqp imp (q). (4.5) 

JQ J-D 

The bulk densities, as was the case for U > 0, are inde- 
pendent of V p . As such the bulk densities can be related 
to those for U > via a spin-j particle-hole transforma- 
tion. This will not be exploited here but is discussed in 
some detail in Ref. HE- 

We can conclude on the basis of this description of 
the U < ground state alone that there are no non- 
analyticities as a function of V p in the various associated 
impurity quantities. V p appears solely in A(g) and A(q) 
is a smooth function of V p . Nonetheless we will compute 
the scattering phases of electrons at the Fermi surface 
and show they have a similar smooth dependence upon 
V p . This will also allow a comparison between our and 
the Dyson equation approach of Section V. 

To determine the scattering phases of electrons off the 
impurity, it is necessary to isolate the impurity portion of 
the momentum. The total momentum is available from 
taking logarithms of the Bethe ansatz equations (|2.33(l 
(up to possible constant terms) 



p(q) 



8=1 



= q+ 6 -^+ f a(X)6i(g(q) — A); 

L JQ 



IQ 

M 



P (x) = 2 -?± = -W{e 2 {x-}!) + 2<K) 
0=1 

i N 

+-^(^(A- s ( ?a )) + 2i) 



<j(X')(6 2 (X-X') + 2ir) 



-D 



dqp(q)(9i(X- g(q))+ 2vr). (4.6) 



The impurity momenta are then found by isolating the 
1/L pieces: 

f° 

Pimp(q) = S(q)+ / (T imp (X)6i(g(q) - X); 
JQ 
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Pimp(A) =-/ a imp (X')(e 2 (X- A') + 2tt) 



+ / dqp imp (q)(9 1 (\~g(q)) + 2ir). (4.7) 

J-D 

We have chosen the branch cuts of the arctans here in 
order to ensure that the impurity momenta take on their 
bare values at the band edges, i.e. Pi mp (<? = ~ oo) = 
—2 tan~ 1 (V p ) and Pi mp (A = — oo) = 0. 

We again have a relationship between the impurity mo- 
menta and the impurity density of states: 

1 



^imp(A) = — dxPimpW; 

Z7T 



2tt 



dqPimp(q)- 



(4.8) 



These relations will allow us to express again the scatter- 
ing phase of an electron off the impurity in terms of the 
impurity occupancy so verifying a variant of the Fricdcl 
sum rule holds for U < and V p finite. 

To compute the scattering phase of an electron one 
must specify the gluing rules between the charge, q, and 
spin, A, excitations. For simplicity let us focus on the 
scattering phase of a spin | electron. In adding a spin | 
electron to the system, we increase the particle number, 
N, by 1 thus forcing an additional real q to be added to 
the ground state. We, in addition, open up a hole in the 
A-distribution. A spin-| electron is then composed of a 
charge g-excitation and a spin A-hole. 

The momentum of the added electron is then a sum 
of the momentum of its component pieces, p e \ = p(k) — 
p(A), and so the corresponding scattering phase (i.e. the 
impurity portion of p e {) is 

6clT =Pimp(&) -Pimp(A). (4.9) 

By virtue of the relations 1)4. 8|l . this can be recast as 



8 e t 



-D 



dqd q p imp + Pimp (A 



-D) 



dQPimp(X) 



2ir dqp imp (q) - 2ir / dAcr imp (A) 
J-D Jq 

-2tan _1 (V p ) 



27rn impT - 2 tan 1 (V P ). 



(4.10) 



Thus we have established a relationship between the scat- 
tering phase and the number of electrons displaced by the 
impurity. 

In the case U > 0, we argued that a discontinuous 
change in the impurity occupancy was mirrored in a dis- 
continuous change in the number of electrons occupy- 
ing the ground state as revealed by exact diagonaliza- 
tion computations. We also argued that the discontinu- 
ity in nj mp as a function of V p was intimately related to 
a U > ground state being composed of two-particle 
bound states. Given that the U < ground state is com- 
posed of a sea of single particle states, how is the lack of a 
discontinuous change in ni mp (U < 0) to be understood? 



Regardless of the sign of U, discontinuous changes in 
the ground state occupancy, n™, as a function of V p 
are seen in the exact diagonalization studies. With the 
U < ground state being composed of a sea of sin- 
gle particle states, changes in n™ are able to be re- 
flected solely in = n bu ik, leaving nf^ p = n imp 
untouched. If An™ = ±1 then An^ k can by ad- 
justed by 1 as the composition of the ground state per- 
mits such changes, in contrast to the U > ground 
state composed of two-particle bound states only per- 
mitting changes satisfying An{^ k = ±2. In this latter 
case we needed to alter An^ p by =Fl in order to satisfy 



An BA 



bulk 



An BA 



An™ = ±1. 



V. CONDUCTANCE FROM A DYSON 
EQUATION ANALYSIS 

In this section we consider an alternative approach 
relying on Dyson equations to compute the linear re- 
sponse conductance. This approach leads to results for 
the conductance, G, considerably different from those of 
the Bet he ansatz presented in this article. We will show 
that this approach, used in Refs. |33 and l3ll is equiva- 
lent to employing the transformation in l|1.16|l by which 
the potential scattering in the Hamiltonian is gauged 
away. We then argue that because the gauge transfor- 
mation incorrectly computes certain finite U quantities, 
so must the Dyson equations. We will further argue that 
the methodology assumes implicitly both V p and U are 
perturbative quantities and so cannot possibly capture 
the non-perturbative physics we believe present in this 
problem. 



A. Overview of Approach 

We will begin by giving a review of the methodology 
used in Refs. |3(] and OH] as few details appeared in the 
two letters. In this approach the current is written as 



1 



-edAriR - n L ) 



i^m,n R }-[H,n L ]}, 



(5.1) 



where n R/L =J2 a J dxc\ R/ L (x)c aR / L {x) and the Hamil- 
tonian, H, is given by l|2.1|l . For the purpose of simplicity, 
we take Vl =V R = V and Vll = Vrr = 0. Thus I can 
be written as 



I = 



V LR J2 J du>(G R +(u;)- 01+^,(5.2) 
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where G + are the standard Keldysh correlators, for ex- 
ample, 

Gm(t) = i(4(<r = 0,i = 0)d(t)); 

G R l(t) = i(c R {x = Q,t = 0)c L {x = 0,t)). (5.3) 

The remaining (RL) Keldysh correlators are given by (so 
as to fix conventions), 

G ffl W = -i6{t)(c L {0,t)c R {0,0)) 

+i0(-t)<4(O,O) Ci (O,t)); 



-ifl(-t)(c i (0 > t)4(0 > 0)) 



+^(i)(4(0,0) Ci (0,i)); 
GfeW = -i(c i (0 J t)4(0,0)). 



(5.4) 



As here we are interested in the zero field transport, we 
suppress all spin indices appearing in the correlators. 

The next step in the development of this approach is to 
recast all of the correlators appearing in (|5.2[) in terms of 
full interacting correlators involving only dot degrees of 
freedom (d/$) together with non-interacting (i.e. V = 
0) correlators of the lead fermions. To this end we can 
write 



G^{V,u)-G- R +{V,u) = V\G]+(0,w){Gu(V,u>) + G&(V,w))-Grf(y,u){^ 

+V (gi+{0, u>)Gfc(V, d) + G fl +(0, u>)G++(V, d) G^(V, u) (G++(0, u) + G^(0, 



+V 2 G++(V,u) 
-V a G+-(V» 
-V 2 G M + (V^) 



-V 2 G++{V,u) 



{G- L + (0, w) + G R +(0, u)) (G++(0, u) + G++(0, u)) 
(Gi+(0, w) + G R +(0, to)) {G R +(0, w) + G R +(0, w)) 
{Gil (0, w) + GH (0, ")) (G++ (0, w) + G++ (0, «)) 



{R<r+L) 



{Gil (0, w) + G^I (0, W )) (G^+(0, w) + G^+(0, w )) - (i£ <-» L) 



(5.5) 



We now must compute the set of lead electron correlators 
at V = (but finite Vp). 

To compute such correlators, we combine the equations 
of motion, 

d t c L = i[H, c L ] = -d x c L - iV LR c R S(x = 0); 



Here / £ and f R are the Fermi functions in the left and 
right leads. Using these mode expansions together with 
the constraints introduced by the equations of motion, 
we find for the G RR correlators, 



d t c R = -d x c R - iV LR c L 6(x = 0), (5.6) 
together with the mode expansions 

c L/R {x,t) = J ^(c +kL/R e ik ^6{x) 

+C- kL / R e lk{x - t) 6(-x) 

c L / R {x = 0,t) = ^{c L/R {x = + ,t) + c L/R {x = 0~,t)). 

(5.7) 

c + k and c_/j mark modes to the right and left of the 
impurity at x = and are governed by 

( c -kL/R c -k'L/R) = 5(k- k')f L/R {u). 



(1 + K 2 ) 2V p 



{l + V 2 ~f R (oj)-V 2 f L (oj)) ; 



G++(0,^) 



2(1 + V 2 ) 2 
GH{0,u) = 



{2f R {u)-l + V 2 {2f L {u)-l)); 



2(1 + v*y 



{2f R (Lu)-l + V 2 (2f L (cj)-l)), 



(5.8) 
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while wc find for the Grl correlators, 

G B +(0,cu) = G+l{0,u>) = G++(0,c) = G^£(0,u) 



V„ 



(± + V 2 ) 



■(/lH-ZaM). 



(5.9) 



Here V p = Vlr/2. consistent with the notation used in 
developing the Bethe ansatz solution. The correlators 
Gll and Glr can be obtained from the above by making 
the substitution L <-> R. 

There is an ambiguity in the result for the Grl corre- 
lators. The Grl correlators satisfy the Dyson relation 

G&M = -VlrY, UGf L {V p ,w)G^ R {V p = 0,u) 

(5.10) 

with a, b = ±. Consistency with these Dyson relations 
demands that G~^£ and G^J are given not by l|5.9|) . but 
as follows 



Grl (0,w) 



(± + V 2 ) 



± 



2(1 + ^)' 



(5.11) 



This ambiguity is perhaps related to the fact that the 
potential scattering term acts at all energies equally and 
has no natural cutoff. 

Substituting these V = correlators f l|5.8|l . I|5.9|l . and 
l|5.11[l ). we arrive at an expression for the current. /, 
given solely in terms of the retarded correlator of the 
dot, G d e d = G dd - G dd : 



I = j I duT(uj)(f L {uj) - /r(u)); 



4V 2 

T( U ) = 



i + v 2 ) 2 

■— ( 
i + v 2 v 



4K^ 



l-K 2 



(1 + K 2 ) 3 



ReG3j(w) 



8V^ 2 

„ f J lrnGggM. (5.12) 



(i + v; 2 ) 



This agrees with Ref. I^ll By taking = /(w ± 

/i/2), the T = linear response conductance is then given 
by 

2e 2 

G= — T(u = 0). 
h 



Dyson Equation Approach and the Bethe 
Ansatz 



In the non-interacting case (U = 0), we can connect 
the expression for the current given in (|5.12(l with the 
scattering phase approach adopted in this paper. At U — 

°> G^iM is s ivcri b y 

D{u) = {{l + V 2 ){u;-e d ) 2 + 2V 2 (LJ-e d )V p + V 4 ); 



And so the conductance takes the form 
e 2 1 1 



(5.13) 



G 



h l + V p 2 D(lu = 0) 

x(l/ 4 + 4t/ 2 e 2 -4^ 2 e d ). (5.14) 

This, as is easily checked, equals the U = conductance 
computed on the basis of scattering phases 



G 



2^sin 2 (i(<5 e -<5 )) 



2— sin' 
h 



^tan" 1 (~~ + V p)+ tan-^Vp)^ . 

(5.15) 



We can also check that a direct computation of the U = 
correlators in (|5.2J) (bypassing the Dyson relations in 
l[5.5[) ) yields an expression for the current consistent with 
the above. Wc have 



1 



1 



I -2V p (Lj-6 d ) 2 -V 2 (u-e d ) 



D(u) 1 + V 2 
\{GR + d ^)-G- d +{u)-{L^R)) = 

-^)TtW( f3+2F ^-^ f )- 



(5.16) 



Substituting the above into (|5.2|l . we again find the same 
expression for the conductance l|5.14[) . 

Thus the two approaches can be made to agree at 
U = 0. Moreover we see that at U = computing corre- 
lators via the Dyson equations is consistent with a direct 
computation (provided we insist on using the V = cor- 
relators found in Q5.11|l ). But what of U ^ 0. To answer 
this question, we recast the two predicted forms of the 
conductance. If we identify 



e = ij(e d + RcI]( W = 0)); 



Y 1 



V 



i + v 2 - 



(5.17) 



where S(w) is the self energy appearing in the full dot 
correlator, 



1 



W - 6 d - E(w) 

1 



w - e d - RcE(w) + iV 2 ' 



(5.18) 



we can recast the expression for the linear response con- 
ductance, G Uyson , arising from the Dyson equations, as 
follows. 



^fDysa 



e 2 W 2 (e + q) 2 
hl + V 2 e 2 + 1 



35 



q 



-—(l - v 2 ) 



(5.19) 



We point out that e, depending only upon the self energy 
of the dot correlator, is actually a quantity computable 
using the Bethe ansatz for the V p = Anderson model. 
The effect of V p upon the dot correlator is to renormalizc 
the bare values of T and e d : 



,-cff 



V 2 



= e d -V p T 



oil' 



(5.20) 



Once this renormalization is taken into account, V p can 
henceforth be ignored. Invoking the Friedcl sum rule, e 
can be expressed as 

e = cot(|n d (t/,ef,l/ off ,y p -0)), (5.21) 

where n d is the number of electrons sitting on a dot with 
effective parameters, e d s and V eS . 

We can recast the conductance arising from the Bethe 
ansatz calculation, 



G = 2^-sin 2 (i( ( 5 e - ( 5 )), 



(5.22) 



in a similar form as l|5.19|l with the identification 



G then becomes 



-cot(2tan _1 (Vy). 

e 2 W 2 (g + q) 2 
hl + Vp 2 e 2 + 1 ' 



G = 2 



(5.23) 



(5.24) 



Given 8 e = 7mi mp — 2tan~ 1 (VJ ) ) and <5 Q = 2tan _1 (VJ,), 
we see that 5 e + 5 = 7mi mp where again ni mp is the total 
number of electrons displaced by the impurity. 
For the two approaches to agree we require 



or 



^imp — 7fil d . 



(5.25) 



Thus we ask that the number of electrons displaced by 
the impurity equal the number of electrons sitting on the 
dot. As we have argued in both Section II and III, this is 
clearly not the case with the Bethe ansatz computation. 
The number of electrons on the dot is a well-behaved 
quantity of both V p and e d whereas rii mp sees numerous 
discontinuities. 

The relations in (|5.25|) make it possible to show that 
the Dyson equation approach is equivalent to gauging 
away the potential scattering term as is done in i|1.16fl . 



With potential scattering gone, the transformed scatter- 
ing phase, <5 e , is simply equal to TTn d by the Friedel sum 
rule. To find the scattering phase of the original electron, 
we undo the gauge transformation giving us 



5 e - 2tan" 1 (V p ) 
Tm d — 2 tan - 1 ( V p ) . 



(5.26) 



Comparing (|5.25|) and l|5.26|l we see by identifying 6 e + S Q 
with irn d , the Dyson equations are implicitly carrying out 
this gauge transformation. 

We now argue that this gauge transformation incor- 
rectly computes the scattering phase at finite U . (On the 
other hand, if one was interested in quantities involving 
the dot degrees of freedom alone, this gauge transfor- 
mation is more than adequate.) We are able to demon- 
strate this at the level of quantum mechanics. Consider 
the two-particle eigenstate constructed in Section II. B. 2. 
This eigenstate is composed of two electrons of energy q 
and p. The scattering phase, 5 q , the electron with energy 
q as it travels from x ~ — oo to x — oo is given by 



Sq(p,e d ,V p ,T) = 



1 (j)(x > 0,q,p,e d ,r) 
i 4>{x < 0, q,p, e d ,T) 
+ S(q,e d ,T,V p ), 



(5.27) 



where <j>{x) is defined in (|2.19|) and H2.20(l while 8{q) is 
given (|2.13|) . The first term in the above involving 4>(x) 
is the phase due to the interaction between the two elec- 
trons. It disappears if U = 0. If the gauge transformation 
led to equivalent results we would find that the scattering 
phase can be written as 

S qiq ,e d ,V p ,T) = -lo %(x<0 ^ £fireff) 



S(q,e 



Off "peff T 7 

d i 1 ! V P 



0) - 2tan" 1 (Vp), (5.28) 



where the final contribution to the phase, — 2tan _1 (Vp), 
is added in undoing the transformation. As is easily 
checked, the two expressions are not equivalent. In the 
presence of finite U, the gauge transformation thus does 
not lead to the correct form of the scattering phase. 

While this is our most definite critique of the Dyson 
equations, we can draw attention to two other more 
amorphous difficulties with the approach. The first stems 
from discrepancies between the direct computation of a 
correlator and computing it using the Dyson equation. 
We have already seen this issue arise at V = in com- 
puting G — (V = 0,cj). And we find it persists to finite 
dot-lead coupling, V. To this end, consider the Dyson 

ret / 
RL > 



relation for Gp^V, w) 



Gg* (0, w) + V 2 (Gf R (0, w) + Gf L (0, w)) 



xG^iV, co) (GT L (0, u) + GT R (0, u)) . 



(5.29) 



If we then employ the correlators in l|5.11|l to determine 
the V = retarded lead-lead correlators appearing in 



36 



we find 



V 2 



1 V P 

2 (1 + V 2 ) 
1 



4(1 + V 2 ) D{uj) 



((uj-e d )(i + V p ) 2 + V 2 (i + V p )). 

(5.30) 



But if we compute G x RL {y,uS) directly using the equa- 
tions of motion l(5.6|) and mode expansions 1)5.7(1 . we find 
instead 



x((, 



V 2 



1 



A{1 + V 2 )D{uj) 
e d )(i + V p ) 2 + V 2 (i 



V p )). (5.31) 



We see that even at V ^ 0, the two methods of computing 
G r R L {uj) give a discrepancy of —V p /2(1 + V 2 ). Agreement 
between the two computations moves much further apart 
if we instead take G^(0, ui) to vanish identically as was 
originally indicated by the equations of motion, i.e. as in 

Thus wc generally find discrepancies between the com- 
putation of correlators using the Dyson relations and 
a direct computation using the equations of motion. 
We could eliminate the discrepancy between (|5.3l)|) and 
(|5.31() by instead (again) computing G RR ^ LL (V, ui) and 

then using the Dyson relations in l|5.10|l to compute G* RL , 
but this is not overly satisfying. The problematic nature 
of the Dyson's relations at U = does not instill confi- 
dence in their use at finite U. 

The second difficulty arises in that the Dyson rela- 
tions express all quantities in terms of the dot correla- 
tor, GJjy (V, w) , and non- interacting (V = 0) lead- lead 
correlators. Now we believe that it is correct to treat 
G t Ij(V,u) as perturbative in both U and V p (just as the 
gauge transformation in l|1.16fl leads to correct results for 
dot quantities). At U = 0, the only effect of V p on G^ 
is to rcnormalize both V and e d - Thus the perturbative 
series in U for the G r d e ^ correlator must be well-behaved 
and convergent for the same reason it is at arbitrary V 
and e d but V p = 0. But while G T ?^{uS) is perturbative in 
U and V p , G r RL need not be (and must not be for consis- 
tency with our solution). Thus expressing G r RL in terms 
of U, V = correlators and G t J^(V,uj) is a mistake. It 
implicitly assumes the problem is perturbative in both 
V p and U. We note that this critique applies only to the 
V p ^ case. If V„ = 0, the reduced Dyson relations (see 
for example Ref. l26t) agree with the exact solution. 

If correct, the origin of this difficulty lies in the manner 
the Dyson relations are established. The Dyson relations 
are not non-perturbative but rather are valid order by 
order in perturbation theory, representing a regrouping 
of terms. The Dyson relations needed at V p ^ (i.e. the 
second relation in (|5.5|) ) involve a regrouping of terms 
at 0{V 2 ) whereas the Dyson relations needed at V p = 
involve a regrouping of terms at only 0(V). It would then 
seem that in this more extensive regrouping of terms that 
non-perturbative information is lost. 
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FIG. 22: U < plots of Trrid, with rid the number of elec- 
trons sitting on the dot, and with the number of 
electrons displaced by the impurity, as a function of ej, — /i, 
where /j, is the chemical potential. In these plots we have used 
U = -1 and r = 0.1. 



C. Dyson Equations and the U < Bethe Ansatz 



While it is obvious the U > Bethe ansatz and the 
Dyson equations give diverging answers for the conduc- 
tance G, it is not a priori obvious the same is true at 
U < 0. At U < 0, as detailed in Section IV, the conduc- 
tance is a smooth function of V p and e^. The quantum 
phase transitions that appear at U > vanish at U < 0. 
Might in this case then the two approaches agree? 

We plot the two quantities, (S e + S a ) = 7rni mp and nna, 
in Figure 22 as a function of t d — /i. We see that at 
U < the number of electrons displaced by the impurity 
is still different than the number of electrons occupying 
the dot. This suggests that recasting the correlators ap- 
pearing in l|5.2|) in terms of the full dot correlators and 
V, U = lead-lead correlators remains inappropriate. As 
our demonstration in Section V.B of differing scattering 
phases between the two approaches was made without 
explicit recourse to the sign of U, this is not surprising. 

We can draw one additional conclusion from the U < 
analysis. The lack of discontinuities/non-analyticities in 
the transport for U < suggests that each term of a 
perturbative expansion of the current in U is well defined, 
but the series itself is not absolutely convergent. If this 
is correct, the Dyson equations are then attempting, in 
some sense, to manipulate the terms in a poorly behaved 
series. 
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VI. COMPARISON WITH GORES ET AL. 

In this section wc study in some detail the experimen- 
tal observations of Fano resonances in the linear response 
conductance through a quantum dot reported by Gores 
ct al. (Ref. ITih . We will focus upon three aspects of 
their observations: i) the general behaviour of the Fano 
resonances as function of gate voltage; ii) their behaviour 
as a function of the dot-lead coupling; and iii) their be- 
haviour in a magnetic field. Before considering any of 
this in detail, one must treat the question of how appro- 
priate it is to compare the predictions arising from the 
model in (|2.1|) with the observations of Gores et al. 

The presence of asymmetrically shaped peaks in the 
conductance reported in Ref. ^] suggests strongly that 
transport through the region containing the quantum dot 
occurs through (at least) two distinct paths, one energy 
dependent and so resonant, one not. However the na- 
ture of the paths is equivocal. In using the Hamilto- 
nian in i|2.1[) . wc model the resonant path as a tunnel- 
ing through the quantum dot while we model the second 
tunneling path in the most minimal way possible: as a 
contact term. The appropriateness of this is not entirely 
clear. The comparison between theory and the experi- 
mental setup would be more robust if the geometry of 
scattering in the dot region was more clear, for example 
if the observed Fano resonances arose from a dot embed- 
ded in an Aharonov-Bohm ring. Notwithstanding this 
uncertainty, we will find that a reasonable comparison 
between our predictions and the observations of Ref. 
can be made. This perhaps suggests that the physics 
is determined merely by the presence of two interfering 
paths and not so much by the exact details of the paths. 

In (|2.1[) we take the dot to have a single spin degenerate 
level. This requires the level broadening, T ~ V 2 , to 
be considerably less than U + Ae, where Ae is the level 
spacing. At least for a subset of the data in Ref. 0, 
this condition is met with T/(U + Ae) ~ 1/35. We note 
that in experimental measurements on dots with a single 
tunneling pathSS, we have T/(U + Ae) ~ 1/12. Here 
the Anderson model does an excellent job of describing 
the scaling behaviour of the reported finite temperature 
linear response conductance. 

An alternative approach to understanding the obser- 
vations of Ref. [Ij was adopted in Clerk et al. (Ref. 



l29f> . There a combination of random matrix theory to- 
gether with a non-interacting Green's function treatment 
modeling the size and shape of the nano-junction (follow- 
ing Ref. lil) was used to describe scattering through a 
dot coupled to a non-resonant channel. Aspects of the 
observations in Ref. ^] were well described by these 
methods. In particular Ref. |2^ was able to reproduce 
generic features of the linear response conductance both 
in a magnetic field as well as its response to changes to 
the dot-lead coupling, T. However these approaches were 
unable to fully explore the consequences of Coulomb in- 
teractions. We, on the other hand, have argued that 
Coulomb interactions lead to highly non-trivial physics 



non-perturbatively related to the non-interacting case. 
Nonetheless it is unclear whether this approach here is 
more appropriate than that adopted in Ref. |2^. At 
least one can say our approaches share some of the same 
physics. We will comment upon this further in the con- 
text of actual comparison with experimental data. 

We now turn to discussing the aforementioned obser- 
vations of Ref. Il4 



A. Behaviour of Fano resonances with gate voltage 

Gores ct. al. report two well developed Fano reso- 
nances. The resonances arc plotted in Figure 23, appear- 
ing as asymmetric dips. These two dips are (slightly) 
differently shaped with different minima and maxima. 

To model these resonances we must first take into ac- 
count the possibility of an asymmetry between the dot- 
lead couplings, i.e. Vl ^ Vr. If Vlr = 0, the only effect 
of Vl 7^ Vr is to reduce the maximal possible conduc- 
tance-: 



h{Vl 



(6.1) 



However with Vlr finite, an asymmetry in Vl ir produces 
a more complex effect. As we indicated in Section II, 
integrability in the presence of an asymmetry, Vr ^ Vl, 
requires an additional term in the Hamiltonian, 



STL = Vllc\c l + Vrrc r cr) 



x=0 



with 



Vll = -Vrr 



Vlr 
2VlVr 



(vl - V 2 ) 



(6.2) 



(6.3) 



Thus in order to take into account the reduction in 
G max due to Vl ^ Vr we must tolerate the presence of 
Vrr / Vll terms. However we believe that Vrr/Vll will 
only quantitatively change the predicted nature of trans- 
port through the impurity. And a qualitative description 
is perhaps the best for which can be hoped given the 
overall simplicity of the model together with the lack of 
complete understanding of the tunneling paths in the ex- 
periment. 

We now consider how to fix the various parameters in 
the Hamiltonian. We take the full width at half maxi- 
mum of the resonances to be indicative of 2T ~ (V 2 +V R ! ). 
As such the ratio of the distance between the two reso- 
nances and r is roughly 1/20. To determine then U/T we 
note that the resonances occur roughly at a dot chemical 
potential corresponding to e™* and e™*. Thus we must 
have 



r 



~crit 



20. 



(6.4) 



Using (|3~T2T) and (|3~53T> we then find that U/T ~ 20 (and 
consequently the distance between the two resonances 
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FIG. 23: A successive pair of Fano resonances observed in Ref. El The black circles mark the experimental data, the solid 
curve the computation from the Bethe ansatz. For this computation, we take U/T = 1/20, V p = 0.78, and Vr/Vl = 0.53. 



corresponds to U). This leaves us to determine Vl/Vr 
and V p . We know from the analysis of Section II that in 
general (unless V p = ±V^ crlt ), the conductance at e d = 
—U/2 and at e d = oo is 



G(e d = -U/2) = G(e d = oo) 



2—7— 
h \l 



W 2 



(6.5) 



For one of the resonances, we then determine the vari- 
ation in the conductance, AG = G ros . max — G rcs . m ; n 
as a function of V p numerically. Knowing AG(V P ) and 
G(e d = oo) as a function of V p permits V p to be then 
fixed. We so find V p ~ 0.78. Alternatively, and more 
crudely, we could use the estimate for AG given in Ij3.16|l . 
If we do so we end up with an estimate of V p = 0.64. The 
difference between the two computations of V p arises in 
that (|3.16|l is only valid for V p large. With V p in hand 
we can then fix the last unknown, the ratio of Vl to 
Vr. From (|6.5|) we find 7 ~ .66 which then leads to 
Vl/Vr ~ .53. 

A final consideration before a comparison can be made 
is the determination of how the gate voltage, V g , applied 
to the dot is related to the dot chemical potential, e d . 



For simplicity we assume V g cx e d . Midway between the 
two resonance marks the point e d = —U/2. We can then 
fix the proportionality constant between V g and e d by 
insisting the positions of the minima of the rightmost 
resonance match. 

Having done this, we are thus able to reproduce the 
behaviour of both Fano resonances in a reasonable fash- 
ion. We do note that the predicted conductance in the 
immediate vicinity of the rightmost resonance overshoots 
the measured conductance by 15%. The predicted con- 
ductance, G, also possesses a number of sharp features 
not seen in the data. However at the particular gate 
voltages at which these occur, we are comparing a zero 
temperature computation with finite but low tempera- 
ture (T = lOOmK) data. As such some of the sharp 
features could be seen as a T = artifact. 

A notable characteristic of this data set is that the 
resonances do not vanish. One possibility is that there is 
some incoherent background to the conductance, some- 
thing assumed in Ref. Il4l An incoherent background is 
necessary if one assumes a non-interacting problem. In 
such a case we may write for the conductance (without 
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function of ed, does not force the vanishing of G. 

Gores et al. report on the temperature dependence of 
the resonances in Ref. They observe that the re- 

gion in between the resonances is markedly insensitive 
to temperature while the resonances themselves grow in 
a Kondo log fashion as the temperature is lowered from 
T = lOOOmK to T = lOOmK. While we are not prepared 
to go into an in-depth treatment of the temperature de- 
pendence here, we will make a few sketchy comments. 

In our discussion of the magnetoconductance in Sec- 
tion III we identified the ranges of values of V p at = 
— U/2 where there would be low- lying spectral weight 
present. Our purpose there was to identify regimes of 
V p where small magnetic fields would drastically change 
the conductance. But what is true about a finite mag- 
netic field is also true of finite temperature. At V v = 0.78 
and with the parameters, U and T, identified as above, 
we expect on the basis of Section III a surfeit of low ly- 
ing spectral weight. Thus as temperature is lowered we 
expect the conductance to vary rapidly in the region of 
€d = —U /2. In Gores ct al. no such variation is observed. 
The cause of this is an extremely low Kondo temperature 
at e d = -U/2. Using2£^ 



FIG. 24: A set of Fano resonances for differing values of V. 
These curves are computed using V p = 0.4 and with Vl = Vr. 



a background contribution) the following, 

e 2 (e + q) 2 



G 



h e 2 + l 



(6.6) 



where e = cot(7rrid/2) (n d is the total number of electrons 
sitting on the dot) and q = — cot(2tan _1 (V / ^)). Thus as 
ed is varied (so changing rid from to 2) , G will necessar- 
ily vanish at some ed- In the Dyson equation approach 
(see Section V) at finite U, G still vanishes for the same 
reason. In either case a background contribution to the 
conductance is needed in order to be consistent with the 
data set. 

As an alternative to explicitly adding a background 
contribution, one can break time reversal symmetry. 
With this symmetry broken, q in H6.6(l is potentially com- 
plex and G then need not vanish as ed is varied. This 
is seen in the random matrix treatment of Ref. |2^ as 
well as in Refs. and where the introduction of 
an Aharonov-Bohm phase serves to break time reversal 
symmetry. 

The approach used in this article, on the other hand, 
does not necessitate either the breaking of time reversal 
symmetry nor the explicit introduction of a Gbackground- 
Rather G remains finite through the non-perturbativc 
effects of U. Although G under the Bethe ansatz can 
be cast in the form (|6.6|l . e in this case does not equal 
cot(7rn^/2). Instead e = cot(7rrij m p/2) where rii mp is the 
total number of electrons displaced by the impurity. This 
quantity, behaving considerably differently than n d as a 



Tk = J EL c n(e d ( td +u)-r 2 )/(2ru) ^ 



(6.7) 



as the expression for the Kondo temperature, we find 

T k {e d = -U/2) ~ lmK. 

As such the temperature range studied in Ref. ^| (100 — 
lOOOmK) is decidedly in the high temperature regime (in 
the region of ed ~ —U/2). Our zero temperature com- 
putation nonetheless finds reasonable agreement because 
the low temperature behaviour, T -C T/., is equivalent 
to the high temperature behaviour, T 3> IV This dif- 
fers markedly from the Anderson model at V p — 0. At 
V p = .78 the low lying spectral weight (as encoded in the 
Pimp (q) distribution) doubles. Crudely speaking then we 
expect the high temperature scattering phase, 5 e , to be- 
have as 

5 e (T » T fe , e d = -U/2) ~ S e (T = 0,e d = -U/2) + 2tt. 

(6.8) 

The shift by 2ir (as opposed to merely it at V p = 0) thus 
implies a coincidence of the high and low temperature 
conductance. It also explains why our T = computa- 
tion adequately conforms to the (high temperature) ob- 
servations (about e d ~ —U/2). 

While at e d = — U/2 the measurements in Ref. are 
high temperature, the corresponding observations in the 
region of the resonances may be considered low temper- 
ature. In the region of the rightmost resonance, 6^ = 
and Xfc ~ IK. This is consistent with the observations 
in Ref. ^] of the resonance growing in depth as tempera- 
ture is lowered. Thus at the lowest temperature at which 
measurements were taken, T = lOOmK, the system is in 
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FIG. 25: The response of a bipolar Fano resonance to the application of small magnetic fields, H . To compute these curves we 
employ T/U = 1/30 (with Vl = Vr) and V p = —3.6. Inset: Observed response in Ref. ^| of a Fano resonance to small applied 
fields. 



a low temperature regime. And consequently we do not 
the behaviour in the two regimes to coincide to obtain a 
reasonable match between data and theory. 



B. Dependence of width of Fano resonances upon T 

In Ref. Fano resonances were studied as a function 
of the total dot-lead coupling strength, T, where it was 
observed that the width of Fano resonances exhibit a non- 
monotonic dependence upon T. (In a dot with a single 
tunneling path, the width of a resonance merely increases 
with r.) Together with this non-monotonicity, the overall 
shape and amount of asymmetry in the Fano resonances 
was observed to be sensitive to the strength of T. 

We can reproduce this array of behaviour. Plotted in 
Figure 24 is the linear response conductance for a set of 
differing Vs. For T small, a Fano resonance appears as 
a sharply peaked bipolar structure. As T is increased, as 
naively expected, the bipolar peak broadens. However at 
some critical value of T ~ (.3 — A)U, the bipolar reso- 
nance is replaced by a narrow unipolar one. With further 
increases in T, this resonance broadens out. 

We point out that the reported behaviour of the Fano 



resonances as a function of T mandates against the form 

S e (q) = -2tan -1 (^— ) - 2tan~ 1 (V;) ! 
q - td 

of the single particle scattering phase (sec l|1.14|l ). This 
is to say, these results imply the scattering due to the 
two paths should be treated as coincident (as contained 
in l|1.12p V This is implicit in both the behaviour of the 
Fano resonances as a function of gate voltage as well as 
the resonances' behaviour in a magnetic field in the fol- 
lowing section. But here we perhaps have a more ex- 
plicit demonstration. If the above was appropriate, vari- 
ations in r would not affect the asymmetry parameter, 
q = — cot(2 tan _1 (V^,)). But as a measure of the shape 
of the resonances, q must and does vary (see Figure 3 of 
Ref. ^3 f° r explicit numbers). 

In Ref. [2^ a different approach is taken to describe this 
phenomenon. Using a specific microscopic picture of the 
dot, these authors consider the effect of modifying the 
second tunneling path (as opposed to directly changing 
the strength of the dot-lead hopping of the first tunnel- 
ing path). They find that the width of the resonances 
behaves non-monotonically as the strength of an impu- 
rity placed in the second pathway is varied. We might 
mimic this approach by varying V p instead of T. If we 
did so, we would certainly find the width of the Fano res- 
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onances varies non-linearly (see for example Figures 12 
and 15). 

C. Linear response conductance at H ^ 

The behaviour of Fano resonances in magnetic fields 
was also studied in Ref. It was found that the res- 
onances exhibited a marked response to extremely small 
magnetic fields (gfiH/T ~ 10~ 2 ). In particular they 
demonstrated that upon application of H , a small bipo- 
lar Fano resonance was transformed into a much larger 
unipolar structure (see inset to Figure 25). The authors 
of Ref. Q suggested the origin of this behaviour lay in 
the large effective area covered by an electron traversing 
the dot via the resonant channel. This large area, A, 
needs a relatively small magnetic field, H, to produce a 
flux A ■ H significant enough to induce the breakdown of 
coherent backscattcring. 

We are able to offer a different explanation found in 
the Kondo-likc physics present in our treatment of the 
problem. From Section III we know that there exist wide 
ranges of V p where a doubling of the low-lying spectral 
weight at positive energies occurs relative to V v = 0. 
This low-lying spectral weight can be accessed by mag- 
netic fields on the order of the Kondo temperature, Tj., as 
computed in the V p = Anderson model. In such ranges 
of V p , and as discussed in detail in Section III, we expect 
a marked sensitivity to small magnetic fields. 

It is then not surprising we are able to reproduce the 
magnetoconductance measurements of Ref. [3 In Figure 
25 we find that the conductance for H = 0, U/T = 30, 
and V p = —3.6 (with this choice of V p the low-lying 
spectral weight has been enhanced in comparison with 
V p = 0) is characterized by a Fano resonance with a 



small bipolar structure. Upon introducing a small field, 
g^H = 0.006r = 1.12Tfc, a unipolar peak is superim- 
posed over the bipolar structure. As the field is increased 
further (although still kept to be the same order as the 
Kondo temperature) this unipolar peak narrows in width. 
Plotted in the inset to Figure 25 are the experimental 
measurements of the conductance in the presence of a 
small magnetic field. At H = 0, Ref. [lj sees a small 
bipolar peak which as a small magnetic field is turned 
on is obliterated by a much larger unipolar structure. 
That the bipolar structure does not entirely disappear in 
our calculation is perhaps a reflection that our computa- 
tion is done at T = 0. With finite temperature blurring, 
the bipolar peak might completely merge into the larger 
unipolar structure. 

In Ref. |2^ this sensitivity to field is explained in a sta- 
tistical sense. The introduction of a magnetic field (of 
any size) serves to break time reversal invariancc. The 
breaking of this discrete symmetry changes the statisti- 
cal S-matrix (i.e. it places it in a different universality 
class) used by Ref. |23 to describe scattering through the 
dot. At least in an average sense this change tends to 
favour larger unipolar over smaller bipolar structures in 
the conductance. 
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